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ABSTRACT 

We study the stellar populations of 1,923 elliptical galaxies at 2 < 0.05 selected from the Sloan 
Digital Sky Survey as a function of velocity dispersion, tr, and environment. Our sample constitutes 
among the largest high-fidelity samples of elliptical galaxies with uniform imaging and optical spec- 
troscopy assembled to date. Confirming previous studies, we find that elliptical galaxies dominate 
at high luminosities (> L*), and that the highest-cr ellipticals favor high-density environments. We 
construct average, high signal-to-noise spectra in bins of a and environment and find the following: 

(1) lower-CT galaxies have a bluer optical continuum and stronger (but still weak) emission lines; (2) 
at fixed a, field ellipticals have a slightly bluer stellar continuum, especially at wavelengths < 4000 
A, and have stronger (but still weak) emission lines compared to their group counterparts, although 
this environmental dependence is strongest for low-cr ellipticals and the highest-cr ellipticals are much 
less affected. Based on Lick indices measured from both the individual and average spectra, wc find 
that: (1) at a given ct, elliptical galaxies in groups have systematically weaker Balmer absorption 
than their field counterparts, although this environmental dependence is most pronounced at low a; 

(2) there is no clear environmental dependence of (Fe), while the a-element absorption indices such 
as Mg b are only slightly stronger in galaxies belonging to rich groups. An analysis based on simple 
stellar populations (SSPs) reveals that more massive elliptical galaxies are older, more metal-rich and 
more strongly a— enhanced. We also find that: (1) the SSP-equivalent ages of galaxies in rich groups 
are, on average, ~ 1 Gyr older than in the field, although once again this effect is strongest at low a; 
(2) galaxies in rich groups have slightly lower [Fe/H] and are marginally more strongly a— enhanced; 
and (3) there is no significant environmental dependence of total metallicity, [Z/H]. Our results are 
generally consistent with stronger low- level recent star formation in field ellipticals at low a, similar 
to recent results based on ultraviolet and infrared observations. We conclude with a brief discussion 
of our results in the context of recent theoretical models of elliptical galaxy formation. 

Subject headings: galaxies: fundamental parameters (classification, colors, luminosities, masses, radii, 
etc.) — galaxies: elliptical and lenticular, cD — galaxies: formation - galaxies: 
evolution - galaxies: stellar content 



1. INTRODUCTION 

The formation and evolution of elliptical galaxies re- 
mains one of the most challenging open problems in 
the general theory of galaxy formation and evolution. 
In the current standard ACDM cosmological model 
(Komatsu et aLl l2010f) . structure grows hierarchically 
(White & Reei [l92i and merging is an unavoidable 
process in galaxy formation. It has long been pro- 
posed that spiral gala xies may even tually merge and 
form elliptical galaxies (|Toomrdll977| ). Recent improve- 
ments in both observations and numerical simulations 
have yielded rema rkable support for the merging picture- 



Deep photometry |F errarese et al."1994; Kormendv et al 
iggi iLauer et al.n'iggS; Kormendv 19991: iLauer et al 



2005; Ferrarcse et al. 2006; Lauer et a l. 2007; Cote et al 



2007tlKormendv et al..,2009l ) with the Hubbl e Space Tele- 



scope (HST) and integ ral-fie ld spectroscopy (iBacon et al.l 
I2001t lEmsellem et al.ll2007t ICappellari et alJi2007D have 



shown that lower-mass elliptical galaxies with M < 
M* in general have cuspy (extra-light) surface bright- 
ness profiles in their centers and are kinematically sup- 
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ported by relatively fast rotation, while more mas- 
sive elliptical galaxies with M > M* have "core- 
like" central surface brightness profiles (i.e., missing 
light), and a.re usually slow rot ators. Recent numer- 
ical studies dM ihos fc Hcrnguist 1994; Cox et al.l 120061 : 
iNaab fc Ostriker 2009 : Hopkins et al..,2009al ) have shown 
that gas-rich mergers between disk galaxies (wet merg- 
ers) can produce a cuspy central surface brightness profile 
and the fast rotational kinematic signature of low-mass 
elliptical galaxies, while subsequent gas-poor mergers be- 
tween less massive elliptical galaxies (dry m ergers) then 
will f orm more massive eh iptical galaxies (jNaab et al.l 
[2006tlHoDkins et al.l[2009hh . 

The mass and luminosity functions of massive red 
galaxies from deep high redshift surveys at z ^ 1 also 
suggest that dry mergers could have played an impor- 
tant role in elliptical galaxy formatio n since redshift one 
(|Ben et al.ll2004bl: iFaber et al.ll2007t ). However, whether 
or not dry mergers are important remains controver- 
sial. Other studies show that massive red galaxies may 
have not undergone many dry mergers since redshift 
unity (jCimatti et al.ll2006t[B7own et al.ll2007l : ICool et all 
I2OO8I ). Merger rate studies also draw vario us conclu- 
sions about t h e significance of dry mergers ([Bell et al 



sions about t h e sigmticance ot dry mergers (laell et al 
[2004al [2006allbl : iMasiedi et ani2006l [2001 iRobaina et al 
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120101) . Meanwhile, studies of early- type galaxies'^ at very 
high redshift show that a significant fraction of massive 
evolved spheroidal stellar systems are already in place at 
very high reds hift (z > 2) Most of them appear to be 
very compact (Daddi ct al. 2005; Longhctti ct al. 200^ 
Toft et all [2007: Truiillo et al. 2007: van Dokkum et al 



20081 ICimatti et al.ll2008l: ISaracco et al.ll2009t ) and it is 
possible that through minor mergers or gradual accretion 
they can evolve to elliptical galaxies at the present day 
(jNaab et al.ll2009l: Ivan Dokkum et al.ll2010D . 

An alternative way to place strong constraints on 
elliptical galaxy formation theory is through the de- 
tailed study of the local universe. Such work has 
been extensively undertaken for decades. Besides 
the deep photometric surface brightness profile and 
integral-field spectroscopic kinematic studies cited above, 
a non-exhaustive list includes: various scaling rela- 
tions (e.g., t he Fabe r- Jackson relation, the fundamental 
plane, etc.) (iFaber fc Jackson KTOTllDressler et al.lll987l : 
Diorgovski & Davis 1987); the color-ma gnitude diagrani 



L973; Bower ct al. 1992; 



20031 



Hogg et al.ll2004l:lBaldrv et al. ll2004iB alogh ct al. lool); 
and a bsorption-l ine in dices (lPeletierll98ikiWorthev et al. | 
1992tlJ0rgensenlll999l: Imager et al.ll20'00at1 Thomas et al.l 



20051 hereafter T05; among many others). These stud- 
ies show that elliptical galaxies have remarkably uniform 
properties and that their stellar mass content is domi- 
nated by old stellar populations. 

The study of environmental effects can also shed light 
on the theory of elliptical galaxy formation and evo- 
lution. It has been known for some time that the 
most m assive and lum i nous galaxies favor high-density 
regio ns (|DressWl980 |: iBinggeh et al.1 119881: iH ogg et al.l 
20031 12004 : Kauffman n et al.l 120041: iCroton e t al. 2005; 



Blanton et al.ii2005at IBlanton fc Mou stakas"2009). Icnd- 



ing support to the hierarchical mod el (ag^, Mq^& White 
[T9M iLemson fc Kauffmannl Il99l iBerlind et al.1 [2001 . 
Numerous studies have also investigated the environ- 
mental dependence of the scaling relations. For ex- 
ample, the fundamental plane for early-type galax- 
ies seems to show a sraall dependence on environ - 
ment fe.g.. lBernardi et al.|[200llLa Barbera et al.ll2010[ ). 
The color-magnitude relations in different environments 
also exhibit a weak, though statistically significa.nt dif- 
ference (e .g., 'Bernardi et al.' '2003b'; 'Hogg etaLj 12004 
ISkibba et a l. 2009; Blanton & Moustakas 200^ A num- 
ber of studies show that the index-velocity disper- 
sion (a) relations such as Mg-cr and Fe-a show no or 
a weak dependence on environment (e.g., J0rgcnse 



se^ 

m 



Il997t iTrager et alll2000at ISanchez-Blazquez et al.ii200 
among others). However, interpretation of these weak 
dependencies is complicated. T05, using a sample of 
124 E/SO galax ies and recent simple stellar populations 
(SSP) models ()Thomas et al.l l2003bl . hereafter TMB), 
find that massive early-type galaxies in low-density en- 
vironments are on average ~ 2 Gyr younger and slightly 
more metal-rich t han their counterpart s in high-density 
environments. Sanchez-Blazquez et al.l ( |2006b ) . using a 
sample of 94 E/SO galaxies reported similar results. 
One surprising recent finding is that many ellipti- 

^ In this paper, we refer to elliptical and lenticular galaxies 
(E/SO) as early-type galaxies, but we only include elliptical galaxies 
(E) in our sample. 



cal galaxies seem to have a non-negligible fraction of 
young stars. Absorption line studies and fundamental 
plane studies favor a "frosting" model in which early- 
type galaxies consist of an old base pop ulation with 
a sma l l amount of younger stars (e.g., iTrager et al.l 
I2000bl : ICMAardt et al.l 120031: iSchiavod 120071 S07 here- 
after). Results from the G alaxy Evolution Explorer 
(GAL EX: iMartin et al.l I2005D . the Spitzer Space Tele- 
scope (|Werner et al.l 120041 ) and the HST have shown 
that a significant fraction of early-type galaxies ex- 
hibit strong ultraviolet (UV) excess, polycyclic aro- 
matic hydrocarbon (PAH) emission and infrared (IR) 
exces s, imply ing possi ble low-level recent star forma- 
tion (lYi et al.|[2 005: Ric h et al.|[2005|:lKavirai et al.ll2007l: 



Schawinski et al. 2007; Temi et al 



, , 120091: lYoung et alj 

2009: Salim fc Rich 2010)7" Recent star formation is 
also consis tent with observed cold gas in many ellipti- 
cals (e.g.. iFaber fc Gallaghei l 11976 1: iKnapp et aD 119851: 
Ivan Gorkom &: Schiminovichr il997l ). ISchawinski et al.) 
(I2007t ) find that the fraction of near-UV (NUV) bright 
early-type galaxies is ~ 25% higher in low-density en- 
vironments, pos sibly due to stronge r low-level recent 
star formation. iKavirai et al.l ()2009t ) suggest that mi- 
nor mergers can account for the inferred amount of 
star f ormation . Usi ng spatially resolved spectroscopy , 
Shapi ro et al.l (|2010l ) (see also Kuntsch ner et all l2010[ ) 
find that star formation in early-type galaxies happens 
exclusively in fast-rotating systems and occurs in two 
distinct modes: the first with widespread young stellar 
populations associated with a high molecular gas con- 
tent , and the second with disk or ring m orphology (see 
also I Youndl2002ll2005l: [Young et al.ll2008l) . They suggest 
the first may be due to minor mergers, and the second 
due to rejuvenation in previously quiescent stellar sys- 
tems. 

The Sloan Digital Sky Survey (SDSS,|^rkiraD[2Q03) 
has provided the largest set of local galaxies with uni- 
form imaging and spectroscopy and offers a great op- 
portunity to study the nearby elliptical galaxies in a 
systematic and homogeneous way. Previous work using 
the SDSS relied on the pipeline outputs in the SDS S 
to select early-type ga laxies fe.g.. iBernardi et al]|2003al : 
lEisenstein et al.l [20031 ). This method does not only se- 
lect elliptical galaxies, but also lenticular galaxies and 
early-type spiral galaxies, which may have different prop- 
erties and formation pathways. Because the SDSS pro- 
vides high-quality imaging for nearby galaxies, we re- 
analyze these images, use their detailed surface bright- 
ness profiles to preselect bulge-dominated galaxies, and 
visually examine the bulge-dominate galaxies to select 
a high-fidelity clean sample of 1,923 elliptical galaxies. 
Our sample constitutes among the largest high-fidelity 
samples of elliptical galaxies with uniform imaging and 
spectroscopy assembled to date. Comparable samples 
are the visually selected earl y-type galaxy sam ples from 
the Galaxy Zoo project (e.g. . ISchawinski et al.| [2009) and 
the Morphologically Select ed Early-types in the SDSS 
proje ct (MOSES, ,Schawinski et al.H2007t iThomas et all 
|201C1() — however, these samples are at larger distances 
than ours (z > 0.05), which at the spatial resolution of 
the SDSS images makes the classification problem more 
difficult. Using this elliptical sample, we study the de- 
pendence on a and environment of the stellar populations 
in elliptical galaxies. 
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The rest of this paper proceeds as follows. In §2, we 
describe our parent sample and our method for selecting 
a high-fidelity sample of elliptical galaxies. In §3, we 
present the photometric properties of our sample and 
determine their local environments. In §4, we study the 
average optical spectra, Lick indices, metallicity and age, 
and their dependence on a and environment. We discuss 
our results in the context of theoretical models in §5 and 
summarize our principal conclusions in §6. 

We adopt a ACDM cosmology with = 0.3, 

17a = 0.7 and i?o = 100 /iioo km s"^ Mpc"^ with 
^100 = 0.7. All apparent magnitudes are on the na- 
tive SDSS photometric system and all absolute magni- 
tudes are on t he AB system and c orrected for Galac- 
tic extinction (Schlcge l et al.l Il998l ) and K-corrections 
(jBlanton fc Roweis.,2007[ ). 

2. DATA AND SAMPLE SELECTION 

2.1. Parent sample 

The SDSS-I and SDSS-II imaged 11,663 square de- 
grees of the sky in ugriz and obtained optical spectra for 
1.6 million objects, ~ 0. 7 million of which are galaxies 



Gunn et al.l 119 98: York et al.l 
Abazaiian et al . 2003). Au- 



with r < 17.77 mag (e.g. , 
[2000t iStrauss et alT[200a 
tomated sof tware performs all of the data processing: 
astrometry (jPier et al.l hOO^ source identifica tion, de- 
blending and photometr y (iLuoton et al. ]| 2001[ ): photo- 
metricity determination (Hogg et al."2001'); calibration 
P^uk ugit a et al.. .1996: Smith et al. 2002); spectroscopic 
targe t selection ( Eiscnstci n et al.l l200it IStrauss et al.l 

§; spectroscopic fiber place- 
; and spectroscopic data re- 



120021: iRichar ds et al., ,200 



2003. 



ment (Blanto n et al.l 
duction. More detail ed descript i ons o f these pipelines 
can be found in Stoug hton et al.l (l2002[ ). An automated 
pipeline called idlspec2d* measures the redshifts and 
classifies the reduced spectra. 

For the purposes of computin g large scale s t ructure 
and galaxy property statistics, iBlanton et al.l (|2005b[ ) 
have assembled a subsample of SDSS galaxies known as 
the NYU Value Added Galaxy Catalog (NYU-VAGC)^ 
We select our parent sample from the NYU-VAGC 
low-redshift catalog, consisting of all galaxies in the 
SDSS with z < 0.05. We use the version of this 
catalog corresponding to SDSS Data Release 6 (DR6, 
lAdeTman-McCarthv et alj l2006t) . which contains 77, 149 
galaxies. 

Due to the difficulty of automatic photometric process- 
ing of big galaxies, the SDSS catalog is missing many 
nearby, bright galaxies, even though they are contained 
within the SDSS imaging footprint; the incompleteness 
begins to become important at r < 14.5. Therefore, 
to ensure a complete parent sample, we include any 
low-redshift galaxies from t he Third Reference Cata- 
log of Bright Galax ies (RC3: lde Vaucouleurs et aLlll991l : 
ICorwin et al.l Il994l ) for which we have ugriz imaging 
from the SDSS, but which are not in the NYU-VAGC. In- 
cluding 10, 474 galaxies from the RC3, our combined par- 
ent sample contains 87, 623 galaxies at z < 0.05 (LowZ; 
see Table p. 

We expect all bona fide elliptical galaxies to occupy the 
red sequence of the color-magnitude diagram, since even 



a modest amount of ongoing star formation will result in 
a blue optical color. Therefore, we isolate galaxies on the 
red sequence using a luminosity-dependent color cut: 

Mg-M,> -0.05 X {Mr + 16.0) + 0.65. (1) 

This cut is very generous in that it extends all the way to 
the edge of the blue cloud, and should therefore include 
most or all of th e blue elliptical galaxies found by re- 
cent studies (e.g. jFukugita et al]|2004l : ISchawinski et al.l 
12001 . There are 37, 026 galaxies that satisfy equation 
([1]), which we define as our photometric red-sequence 
parent sample (PhotoRS). Within this sample, 32,726 
have spectroscopy from the SDSS, which we define as 
our spectroscopic red-sequence parent sample (SpecRS; 
see Table [J). 

2.2. Elliptical sample selection 
Traditionally, galaxies are classified int o different mor- 



195S 


■ 


2007 





http : //spectro . pr Inc eton . edu/idlspec2d_ Install . html | 
http: //sdss .physics . nyu. edu/vagc/ I 



Lintott ct al. 200q). Unfortunately, visual classifi- 
cation is inherently subjective, and not feasible for sam- 
ples consisting of tens of thousands of galaxies, which has 
led to a concerted effort by many different groups to clas- 
sify galaxies us i ng ob j ective, quantitative criteria (e.g., 
Consehce et al.l 120031: iLotz et al.l 12004 IScarlata et ahl 
2007^. However, although quantitative morphological 
classification schemes have become increasingly sophis- 
ticated, they arguably do not capture the detailed varia- 
tion in morphological signatures a pparent to the trained 
eye ( see the recent discussion in IBlanton fc MoustakasI 
I2OO9I and references therein). Moreover, in studies of 
elliptical galaxies, contamination from early-type disk 
galaxies (i.e., SO/Sa), can be significant. Indeed, ob- 
jective classification methods have tremendous difficulty 
distinguishing bulge-dominated disk galaxies from true 
elliptical galaxies, although the two galaxy types may 
have experienced very different formation pathways. 

Given these various issues, we opt for a hybrid ap- 
proach, in which we conservatively preselect elliptical 
galaxies using well-defined quantitative criteria, and then 
use visual inspection to remove contaminants. Thus we 
retain the objectivity of quantitative methods and the 
ability to analyze a large parent sample, while still re- 
lying on the superb ability of the eye to identify subtle 
but important morphological signatures like faint spiral 
arms, dust lanes, bars, etc., to help remove contaminants. 
In the next three sections we describe our elliptical galaxy 
selection in more detail. 

2.2.1. The Ellipse method 

We first reduce the two-dimensional (2D) images to 
one-dimensional (ID) sur face brightne s s pro files using 
the Ellipse algorithm (iJedrzeiewskil Il987l see also 
lYoung et all Il979l : iKend |1983D . The Ellipse method 
is widely used in the imaging analysis of galaxies be- 
cause the isophotes of galaxies are well approximated by 
ellipses. 

The basic idea of the Ellipse method is to expand 
the intensity aroun d an ellipse in a Fourier series (see 
iBender et all Il989l for the expansion in polar coordi- 
nates) : 

I{E) = la + Y^ [An cos{nE) + B„ sm{nE)] , (2) 

n=l 
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where / indicates the intensity, E is the eccentric 
anomaly, and An and Bn are the higher-order Fourier 
coefficients. For a given major axis with length a, we use 
n = 1 and n = 2 to find the ellipse that best matches 
the measured isophote. Specifically, we first provide an 
initial guess for the center of the ellipse (xq, yo), the ellip- 
ticity (e), the position angle {(f)), and the intensity (/q). 
We then calculate An and Bn (n = 1, 2) and update the 
parameters by calculatin g the expected dev iation from 
the true parameters fsee iJedrzeiewskil I1987D . and iter- 
ate until the maximum of An and Bn {n = 1,2) is less 
than 4% of the intensity (/q). Next, we measure the 
third- and fourth-order harmonics of the resulting in- 
tensity distribution using least-squares minimization. If 
an isophote is a perfect ellipse then all the coefficients, 
{An, Bn), n = 1, • • •, oo will be identically zero. Non-zero 
coefficients indicate the amount by which the isophote 
deviates from the shape of an ellipse. In particular, 
A4 > indicates a "disky" isophot e, while A4 < cor re- 
sponds to a "boxy" isophote (e.g.. iBender et ahlllQSgf ). 

We apply the Ellipse method to each galaxy image 
in all five bandpasses as a function of the major axis 
radius, using a one-pixel step size (0.396 arcsec) along 
the major axis. At each step, we save all the parameters, 
including the position of the origin (xo,yo)i the major 
axis length (a), the ellipticity (e), the position angle (0), 
the intensity (/o)j the harmonic coefficients {An and Bn 
with n = 1,2,3,4), and the total flux within each ellipse. 

2.2.2. Bulge-disk decomposition 

Given the Ellipse fitting results, we model the sur- 
face brightness profile, Irj, in each bandpass j — g,r,i 
simultaneously as a function of radius R using a two- 
component "bulge plus disk" model: 

Irj = Ibj exp{-KBR^) + Id.j cxp{-KdR), (3) 

where the first term represents a Sersic model of index 
n fe.g.. lSersid[T96llCiottil[T99llGraham fc Drivedl2005l) 
for the bulge component {B), and the second term is 
an exponential model for the disk component (D). The 
parameters Kb and Kd are length-scale factors for the 
bulge and disk component, respectively; they are con- 
strained to be equal in all three bandpasses, while the 
amplitudes Ibj and lo.j are allowed to vary. Note that 
we have chosen not to include the u- and 2:-band surface 
brightness profiles in this analysis because the photome- 
try in those bands is considerably noisier. 

2.2.3. Preselection and visual inspection 

Our goal in this section is to select a high-fidelity sam- 
ple of elliptical galaxies. As discussed in tj2.21 our strat- 
egy is to apply conservative parameter cuts to our ellipse- 
fitting results to preselect a sample of spheroidal galax- 
ies, and then to visually inspect the resulting sample to 
remove contaminating disk-dominated galaxies. 

First, we require the bulge-to-total ratio, B/{B + D), 
to be larger than 0.7 in all of the three bandpasses. 
We then reject galaxies with ellipticities larger than 0.6, 
which only rejects a handful of real elliptical galaxies 
(see §3.1 and Fig. HJ. Next, we parametrize the "fea- 
turelessness" of the surface brightness profile by fitting 
a straight line to the ellipticity versus radius profile. 



and compute assuming uniform errors for the ellip- 
ticities. We define a featurelessness parameter by mul- 
tiplying the by the largest ellipticity in the pro- 
file. By exploring training sets, we apply a generous 
cut in this parameter to reject non-elliptical contami- 
nants while simultaneously minimizing the number of 
real elliptical galaxies that are excluded. Finally, because 
the instrumental dispersion of the SDSS spectrograph is 
69 km s~^, a velocity dispersion measurement smaller 
than 70 km s~^ is not reliable^; therefore, we only in- 
clude galaxies with cr > 70 km s~^. Among the 32, 726 
galaxies in our SpecRS sample (see Tabled]), 22,621 ob- 
jects have a > 70 km s^^, of which 2,648 survive our 
preselection cuts. 

Figure [1] presents color images of a randomly selected 
subset of our elliptical galaxy sample, sorted by increas- 
ing velocity dispersion. Given our generous parame- 
ter cuts, however, our procedure does result in some 
non-elliptical contaminants, three examples of which are 
shown in Figure [2] Among the contaminants are bulge- 
dominated SBO galaxies with very small bars (e.g.. Fig. [2] 
left), galaxies with faint dust lanes (e.g., Fig.^ middle), 
and SO galaxies with faded spiral structures (e.g.. Fig. [2] 
right). Therefore, we visually examined the preselected 
elliptical galaxy sample and excluded such galaxies. Note 
that in Figure [T] we have already excluded such galaxies. 

Our final sample (Elliptical; see Table [1]) contains 
1,923 elliptical galaxies with tr > 70 km s~ . For those 
without SDSS spectroscopy, we have performed the same 
analysis and selected 430 elliptical galaxies. We refer to 
them as the bonus elliptical sample (Bonus). Though we 
do not include them in the analysis in the rest of this 
paper, we make them publicly available, along with the 
Elliptical sample with the SDSS spectroscopy. 

Our selection certainly still suffers from some subjec- 
tivity and may not be "complete" depending on one's 
definition of an elliptical. Because we aim for a clean 
sample, our visual inspection is strict, and we likely ex- 
clude some elliptical galaxies that appear ambiguous. In 
addition, the bulge and disk decomposition classifica- 
tion also possibly identifies real elliptical galaxies with 
low Sersic index (low concentration) as disk-dominated 
galaxies. 

For comparison, if we selected early-type ga lax- 
ies in the same way as in iBernardi et al.l ()2003al ) or 
lEisenstein et al.l ()2003[ ). we would select ~ 5000 galax- 
ies, including ~ 1600 of the galaxies in our Elliptical 
sample. Therefore, we could select most of galaxies in 
our sample with their selection method, but lenticulars 
and early-type spirals would comprise about two thirds 
of the sample. 

3. THE ELLIPTICAL SAMPLE 

Our final sample contains 1, 923 elliptical galaxies with 
(7 > 70 km s^^ and z < 0.05. In this section, we present 
the general properties of the sample, compare them with 
the parent and red-sequence samples, and define the local 
environment of each galaxy. 

3.1. Nature of the red sequence 

^ |http : //www . sdss ■ org/dr7/algor Ithms/veldlsp ■ html | 
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Galaxies can be broadly divided into two groups, the 
red sequence and the blue cloud, acco rding to their 
broad-band color and luminosity fe.g. iBlanton et alj 
I2003bt iKauffmann eTaD [200l iBaldry et all 120041) . In 
general terms, red galaxies are early-types with a 
spheroidal surface brightness profile, that are dominated 
by an old stellar population. However, disk galaxies that 
are dusty or have star formation quenched can also ap- 
pear red in optical broad-band colors. Thus, the red 
sequence in general consists of both disk-dominated and 
bulge-dominated galaxies. 

For example. Figure [3] shows the broad-band properties 
of the Elliptical sample (red dots) compared with the 
parent samples. The top left panel of Figure [3] shows 
the full LowZ sample, in which elliptical galaxies (red 
dots) clearly dominat e the luminous end L*) of the 
red s equence (see also lMarinoni et aLlll999tlBundy et al.l 
l2009f ). The diagonal line corresponds to equation ([1]), 
the color cut used to define the red sequence in §2.11 

In Figure [3] (top-right) we examine the red-sequence 
galaxies (SpecRS) in more detail by plotting the bulge- 
to-total ratio {B/ (B + D)) versus cr. Red-sequence galax- 
ies naturally separate into disk- and bulge-dominated 
groups. Note that our cuts in a and B/{B + D) are vis- 
ible in the distribution of the Elliptical sample. Also 
note that our cut on bulge-to-total ratio tends to exclude 
dwarf elliptical galaxies, whi ch often exhib i t exponen- 
tial surface brightness profiles (lMichardlll985t [Schombert 
1386; Prugnicl & Simicn 19971 ICaon et al.lll993l : lGraham 
2002c iKormendv et al.ll2009l ). 

Although we do not use it for classification, we also ex- 
amine the color gradients of the galaxies in our sample. 
The lower panels of Figure [3] show the color gradients of 
the red-sequence galaxies as a function of Mr (PhotoRS, 
left) and a (SpecRS, right). We define the color gradi- 
ent as the color difference between the g — i color within 
the 15% light radius and that between the 15% and 90% 
light radii. We find that galaxies form two distinctly 
different sequences. The central regions of most galax- 
ies are redder than their outer regions. For the disk- 
dominate d galaxies the color gradient increases with lu- 
minosity ([Mi^iEin [llM Eaiiniran [200^ Mean- 
while, for the elliptical galaxies, t he color gradient de- 
cre ases marginally w ith luminosity (jBoroson et al.l[l983l : 
see ISuh et ahl |2010| for an extensive study in a sample 
similar to ours). These results suggest that the color gra- 
dient can be us ed as an indicator for morphological clas- 
sification (e.g.. iPark fc Choi.2005<) . The physical cause 
of the color gradient is likely that the metallic i ty de- 
creases w ith increasing radius . iKuntschner et al.l (|2010l ) 
(see also ISpolaor et al.l l2009f ) show that for low-mass, 
fast-rotating early-type galaxies the metallicity gradient 
increases with mass, while for more massive systems the 
metallicity gradient becomes shallower, leading to the 
most massive systems being slow rotators with relatively 
shallow metallicity gradients. This result is consistent 
with the ff/luminosity-color gradient relations we find 
here. 

In the top panels of Figure S] we compare the distri- 
butions of tJ, and ellipticity (e) for the Elliptical and 
SpecRS samples. Because e in general varies within a 
galaxy, we must choose a particular definition; here we 
use the mean value between the 30% and 60% light radii. 
We find the distributions of both quantities to be very 



different among the two samples. Specifically, elliptical 
galaxies contain many more massive galaxies, and have 
a median e ~ 0.2. We also note from Figme^{top-right) 
that the e < 0.6 cut we applied in i j2.2.3l reiects very few 
real elliptical galaxies. 

The bottom panels of Figure |4] quantify the fraction 
of elliptical galaxies in the SpecRS sample as a func- 
tion of a and e. We find that the fraction varies from 
< 10% for a < 100 km s"^ and e > 0.5, to ~ 80% 
for a > 300 km s^^ and e < 0.1. Note that our selec- 
tion was intended to be strict; therefore these fractions 
should be treated as lower limits. Nevertheless, these re- 
sults suggest that only if one restricts to the luminous end 
are red-sequence galaxies typically giant elliptical galax- 
ies. At fainter luminosities (L < L*), disk-dominated 
galaxies constitute the bulk of the red sequence. As we 
pointed out above, part of this trend with luminosity is 
due to our bulge-to-total cut; lower luminosity ellipticals 
in general have smaller Sersic indices (e.g., Graham 2002]; 
iKormendv et al.ll2009f ) and our procedure may start ex- 
cluding those when they become consistent with expo- 
nential profiles. Selecting a reliable sample of fainter el- 
liptical galaxies would require more information than we 
use here. 

3.2. Local environment 

There are various ways to define the loc al environment 
of a galaxy (iBernardi et al1[2"0 03a: Eisenst ein et al.ll2003l: 
iBalogh et all 12004 iPark et al., .2007; .Schawinski et al.l 
2()()7). For instance, parametrizing the eirvironment with 
local density separates the central and the outer regions 
of a cluster, while using group membership can separate 
galaxies associated with groups from those isolated in 
the field. Here, we will compare galaxies in groups to 
isolated field galaxies. 

We use a Friends-of-Friends (FoF) grouping algorithm 
to define groups. We choose linking lengths to be 
0.88 Mpc and 1.76 Mpc (120 km s"^) in the angular and 
line-of-sight directions, respectively. We choose the link- 
ing lengths to be short enough to break filaments into 
knots, but not too short to only select the galaxies in the 
center of any groups. The line-of-sight linking length is 
small com pared to other implementation s in the litera- 
ture (e.g., LG oto"200jt iBerlind et al.ll2006[ ): for example, 
iBerlind et al. (2006) use 300 km s ^. We explain below 
why we make this choice. 

We define the groups based on the LowZ sample, which 
includes the RC3 catalog galaxies necessary to ensure 
completeness at the bright end. However, consistently 
finding groups across redshift requires a volume-limited 
sample. At the faint end, the fiux limit of the SDSS 
spectroscopic survey is r = 17.77, which corresponds 
to Air ^ —19 at z = 0.05. We therefore only include 
galaxies in LowZ brighter than Mr = —19 in our FoF 
analysis. The resulting environmental sample (Environ) 
has 57,885 galaxies in total. In our Elliptical sam- 
ple, there are only 16 galaxies fainter than Mr = —19; 
the velocity dispersion cut thus has provided roughly a 
volume-limited sample. In the following sections, we do 
not consider these 16 faint ellipticals when studying stel- 
lar populations. 

We divide galaxies into three broad categories: rich 
group, poor group, and field. Rich-group galaxies are in 
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groups with at least 5 members, poor-group galaxies are 
in ones with 2 — 4 members, and field galaxies are the sole 
members of their group. We note that the line-of-sight 
linking length (120 km s~^) we adopt is relatively small 
compared to the typical velocity dispersion of clusters. 
This choice may result in identifying a galaxy associated 
with a group as an isolated galaxy. Therefore for each 
rich group we define an ellipsoid with radius 0.88 Mpc in 
the angular direction and radius 360 km s~^ in the line- 
of-sight direction. For each poor group, we also define 
such an ellipsoid but adopting a slightly smaller radius 
240 km s^^ in the line-of-sight direction. If an isolated 
galaxy lies within the ellipsoid of any group, we consider 
its group classification to be ambiguous. There are 129 
ambiguous galaxies in total, 46 of them with possible 
association with rich groups and 83 with poor groups. 
Including them either in the sample of field galaxies or 
group galaxies does not affect any of our results; there- 
fore, we exclude them from further analysis. 

In m we study the environmental dependence of the 
stellar populations of our sample. However, it is im- 
portant to exclude low signal-to-noise ratio {S/N) spec- 
tra and strong AGN. We thus will not consider galax- 
ies with spectra of median S/N (per pixel) less than 
10. To exclude strong A GN, we use a standa rd emission 
line diagnostic diagram ()Baldwin et al.|[T98lL BPT). For 
galaxies with emission lines H/3, [O III] A5007, Ha and 
[N II] A6584 of S/N > 5, we exclude them if they ap- 
pear to the righ t of th e demarcation line defined by 
iKauffmann et al.l ()2003[ ) in the BPT diagram. We ex- 
clude 39 strong AGN in this process, which we have ver- 
ified does not affect the analysis presented below. 

To summarize, after excluding 59 galaxies with Mr > 
-19 (16), spectra of low S/N (4), strong AGN (39), and 
129 galaxies with ambiguous group association, out final 
Elliptical sample consists of 347 field elliptical galax- 
ies, 682 poor-group galaxies, and 706 rich-group galaxies. 

Figure [5] shows the angular distribution of the 
Elliptical sample. In the top panel, we show the distri- 
bution of the red-sequence galaxies in the parent SpecRS 
sample. In the bottom panel, we show the distribution of 
part of the Elliptical sample, where we represent field 
galaxies as blue open triangles, and rich-group galaxies 
as magenta small solid squares. We also show the posi- 
tion of the largest group, the Coma cluster as a large red 
solid square. 

4. STELLAR POPULATIONS OF ELLIPTICAL GALAXIES 

In this section, we study the stellar population of our 
Elliptical sample as a function of a and environment 
with average optical spectra and Lick indices. 

4.1. Average spectra 
4.1.1. Stacking method 

We stack the spectra in various samples of velocity dis- 
persion and group classification, as described through- 
out. There are three versions of stacks we make here: 
unsmoothed, smoothed, and uniform stacks. 

For each sample, we create at least the unsmoothed 
and smoothed average spectra. In the unsmoothed ver- 
sion, we do not smooth individual spectra before aver- 
aging them. In the smoothed version, we first smooth 
each spectrum to the largest <j in that sample. The un- 



smoothed version allows us to see more details that are 
washed out in the smoothed version, while the smoothed 
version allows us to compare samples across different en- 
vironments consistently. 

For samples in different a bins, we also calculate a 
uniform version by smoothing each individual spectrum 
to 325 km s^^, to match the largest a in the whole sam- 
ple. These stacked spectra allows us to compare different 
samples as a function of velocity dispersion consistently. 
We use this uniform version to measure the Lick indices 
for the average spectra in the following sections (§4.2.3 
and §4.2.4). 

Figure [5] shows the velocity dispersion distribution as 
a function of group classification. The galaxies in rich 
groups tend to have higher a than those in low-density 
environments. As we show explicitly below, and is al- 
ready well known, the properties of elliptical galaxies de- 
pend strongly on a and luminosity. For instance, the 
color of elliptical galaxies is redder at hi gher luminosity 
and larger a (Figure H e.g., 'Baum' 1959). 

Thus, to make a direct comparison between samples 
in different environments, we need to ensure that we 
compare samples with the same a distribution — oth- 
erwise we will simply be measuring the dependence on 
a itself, which dwarfs all other effects. To do so, we 
weight the galaxies as a function of a in order to achieve 
the same effective cr distribution for each subsample. 
Within each environmental subsample, for a galaxy with 
a certain a, we determine the number of galaxies within 
Act = ±10 km s~^ in both the full sample and in the 
subsample. We then weight each spectrum in the stack 
by the number in the full sample divided by the number 
in the subsample. In this way, the different environmen- 
tal samples have the same effective a distribution. This 
procedure allows us to separate the environmental de- 
pendence from the a dependence. 

The SDSS spectra are observed between wavelengths of 
~ 3800 A and - 9200 A. Because our redshift limit is 2: = 
0.05, for consistency we calculate the average spectra in 
the rest-frame wavelength range available for all galaxies, 
from 3800 A to 8800 A. Before stacking the spectra, we 
also normalize each spectrum to its mean flux between 
5200 A and 5800 A, where the spectra are relatively fiat. 

The spectrophotometry^ in the SDSS is calibrated with 
F8 subdwarfs and is accurate at the few percent level. In 
any case, our differential results between different envi- 
ronments and velocity dispersion samples are largely un- 
affected by systematic errors, unless their strength cor- 
relates with those parameters, which is unlikely. 

4.1.2. Velocity dtsperston dependence 

Figure [7] shows the average unsmoothed spectra as a 
function of a. We overplot two horizontal lines at 0.95 
and 0.35 in each panel to guide the eye. We define the a 
bins to be roughly equal size in log space. The three bins 
(in logioCT) are [1.84, 2.10], [2.10, 2.30] and [2.30, 2.50]. 
The median cr in each bin are 100, 167 and 232 km s^^, 
respectively. The average spectra at all cr are typical of 
an old stellar population, with strong absorption features 
such as the 4000 A break (mainly caused by the Call H 
and K lines), the G-band at 4300 A, the 5180 A Mgl and 

^ |http: //www. sdss . org/dr6/algorithms/spectrophotometry .html| 
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H feature, the 5890 A Na I and the 7200 A TiO lines. 
Though we have excluded strong AGN when calculat- 
ing the average spectra, weak emission features such as 
Ha are still visible, though at a fairly weak level. They 
can either be caused by low-level star formation or weak 
AGN activity. Line ratio diagnostics show that some of 
the ellipti cal galaxies with emission are LINE R- like (see 
also, e.g. JYan et"aIll2006l:[G7aves et al.ll2007t) while oth- 
ers show signs of ongoing star formation. 

Not surprisingly, we see a clear dependence on ve- 
locity dispersion of the average spectra (which is why 
correcting for the a distribution is so important when 
studying environmental effects). The average spectra in 
lower a bins are relatively bluer than those in higher a 
bins, whic h is consistent with the color-m agnitu de rela- 
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, among many others). We examine this result in 



more detail in Figure El where we show the ratio of the 
average spectra (uniform version) of all galaxies within 
different a bins to that of the whole sample. We find 
that lower-tT galaxies exhibit a bluer continuum, weaker 
metal absorption features and slightly stronger emission 
lines (e.g.. Ha). 

However, the cause of the bluer continuum for elliptical 
galaxies at lower a is less clear. The primary difficulty in 
interpretation is the well-known age-metal licity degen- 
eracy in stellar population a nalys is fe.g.. | Fab eij |1972l, 
19731 Oconncll 1980; Rosj |1985|: iRenzini fc B uzzoni 
T98fit iWorthev et al.l Il994t iBruzual fc GhaTIoti l2003f ). 
Both increasing age and increasing metallicity have ex- 
tremely similar effects on the spectral energy distri- 
butions (SEDs). This similarity means that either a 
younger population or a lower metallicity population 
(or some combination) can be responsible for the bluer 
continua of elliptical galaxies at lower a. Further- 
more, the presen ce of metal-rich, old blue hor i zonta l 
branch stars (e.g., ISweigartlll987l: iLee et al.lfTMol I1994D 
can mimic a younger stellar population and introduce 
more flux at t he blue end as well (e.g., 'Yi et al. l 119971: 
Maraston fc T homas 2000; Lee et al. 2000; Conro v et al. 
2010") Blue stragglers fe.g.. IBailvnl 119951: IBrown et al. 
2005,) that extend blueward and brightward of the main- 
sequence turnoff point may also ha ve a non-negligible 
effect on integrated spectra (e.g., IXin fc Dend l2005t 
IConrov eF al. 2010). Stronger AGN activity or less dust 
extinction can also cause the same effect. 

Because of these effects, we cannot be certain what is 
responsible for the velocity dispersion dependence in the 
average spectra. However, the stronger emission lines, 
part of which certainly is due to star formation activity, 
and the weaker metal absorption features suggest that 
both younger ages and lower metallicities are responsible 
for the bluer continuum at lower a. 

4.1.3. Environmental dependence 

Figure |9] shows the average unsmoothed spectra as a 
function of environment. Recall, in these stacks we have 
weighted the spectra so that the effective a distributions 
are the same; otherwise, stronger trends would be visible 
due to the correlation of a with environment. The aver- 
age spectra for all the environmental samples look strik- 
ingly similar to each other. These spectra reveal that 
old stellar populations dominate the stellar components 



in elliptical galaxies in all environments. 

Although by eye the average spectra for all the sam- 
ples look exactly the same, careful inspection shows sub- 
tle variations. In Figure 1101 we present the ratio of the 
average spectra (smoothed version) of different samples 
in different environments to that of the whole sample. 
We find weak but significant differences between the av- 
erage spectra in different environments. Compared to 
the average spectra of the whole sample, the field galax- 
ies have a bluer continuum, especially at wavelength 
< 4000 A, by ~ 1 percent, and have stronger (but still 
weak) emission lines of H/3, [O III] AA4959, 5007, Ha, 
[N II] AA6548,6583 and [S II] AA6717,6731 etc. The 
rich-group galaxies, on the other hand, have a redder 
continuum and have less flux at the blue end by ^ 1 
percent, and have weaker emission lines. 

To further inspect the environmental dependence, we 
compare galaxies in different environments in each bin of 
velocity dispersion. In Figure 111! we show for each a bin 
the ratio of the average spectrum (uniform version) of the 
field galaxies to that of the rich-group galaxies. Interest- 
ingly, we see that in all three bins, the average spectra 
of the field galaxies are bluer than their rich-group coun- 
terparts. However, the weak environmental dependence 
appears to be a strong function of a. In the lowest a 
bin, the field galaxies have about 5% more fiux at the 
blue end than the rich-group galaxies. In the highest a 
bin, the field galaxies have only about ~ 1 percent more 
flux at the blue end than the rich-group galaxies. The 
stronger emission lines in the low-cr field galaxies also 
appear to vanish for the high-a field galaxies. 

As was the case for the dependence on cr, the cause 
of the weak environmental dependence of the average 
spectra is not absolutely clear. However, the difference in 
emission lines suggests that stronger (but still low-level) 
recent star formation could be the cause of the bluer 
spectra in the field galaxies. In the next subsection, we 
examine this hypothesis more carefully. 

4.1.4. Discussion: stronger low-level recent star formation 
activity at lower velocity dispersion and in the field? 

Physical degeneracies make it difficult to derive ac- 
curate stellar population properties from the spectra of 
galaxies with certainty — such attempts suffer from the 
age-metallicity degeneracy, the possibility of poorly un- 
derstood blue horizontal branch and/or blue straggler 
populations, and other effects. In this section we take a 
simplified view and evaluate what kind of young stellar 
population would be necessary to explain the observed 
trend of SED with environment. 

First, we assume that elliptical galaxies consist of 
an old base stellar population and a small frosting of 
young stars from low-level recent star f ormation (e.g., 
iTrager et al.1l2000aHGebhardt et al.ll2003l: S07). Further- 
more, we assume that the strength of the recent star for- 
mation is the main factor that causes the a dependence 
and the environmental dependence. 

Under these assumptions, we fit each individual spec- 
trum with a model consisting of a simple combination 
of old and young stellar populations. The old compo- 
nents consist of 3 old stellar populations with ages and 
metallicities as follows: 13.5 Gyr and solar abundance 
(Z©); 15 Gyr and 0.4 Z©; and 12 Gyr and 2.5 Zq. The 
young components consist of 3 young templates with 
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ages and metallicities as follows: 0.9 Gyr and Z©, 1 
Gyr and 0.4 Z©, and 0.8 Gy r and 2.5 Z© respe c tively . 
All these templates are from iBruzual fc CharlotI (l200l 
SSP models with the Padova (1994) library of stellar 
evolution tracks (e.g., [A longi ct al.' 1993; Brcssa n et al.l 
1993HFagotto et al.lll994a , h : Girardi et al..,1996 ') and the 
Chabriei](|2003l i" initial mass function from 0.1 to 100 Mq. 
For each individual galaxy spectrum, we mask out the 
emission line regions and fit the continuum to a nonneg- 
ative linear combination of these six templates. 

We evaluate how much recent star formation there is 
by comparing the resulting mass fraction fy contained in 
the young templates in the fits. We choose fy = 2% to be 
the minimum fraction of the young component at which 
we consider the recent star formation to be significant. 
This choice allows a convenient comparison with previous 
work using GALEX NUV photometry. For reference, 
fy = 2% gives the color NUV— r ^ 5.4, w hich is the 
color cut adopted bv lSchawinski et all (|2007[ ) to indicate 
recent star formation. 

In Figure [T^ we show that the fraction of galaxies 
with significant recent star formation {fy > 2%) is a 
strong function of both environment and velocity dis- 
persion. Overall, about 20% of the elliptical galaxies 
have more than 2% of the mass contributed from the 
young components. In the lowest tr bin, about half of 
the elliptical galaxies have more than 2% young compo- 
nents. Meanwhile, in the highest a bin, only < 10% have 
fy > 2%. The fraction of galaxies with recent star forma- 
tion also appears to be larger for the field galaxies than 
the rich-group galaxies, by 30% in the lowest a bin. 
As suggested by the average spectra, the environmental 
difference vanishes in the highest a bin (Figure [TT|) . 

Interestingly, our r esults agree reasonably well with 
previous work (e.g., ISchawinski et al.l l2007f) . Using 
GALEX NUV photometry of a sample of 839 bright 
early-type galaxies with Mr < —21.5, most of which are 
elliptical galaxies, sele cted in the SPSS in th e redshift 
range 0.05 < z < 0.10, ISchawinski et al.l ()2007t) find that 
30% ± 3% of massive early-type galaxies have NUV— r 
color bluer than 5.4. They also show that the fraction 
is 25% higher in the lowest-density environment, which 
is defined based on the number and distance of galaxies 
around each object (the density). Despite the differences 
in the definition of environment, as well as the different 
samples and data used, our results are in fair agreement 
with each other. However, they did not see the luminos- 
ity/mass dependence that we see here, probably because 
of their smaller sample size and luminosity range. 

Using a sample o f early-type ga laxies at 0.05 < z < 
0.06 from MOSES, iThomas et all (pOlO) find that the 
distribution of ages is bimodal with a strong peak at old 
ages and a secondary peak ~ 2.5 Gyr younger containing 
~ 10 percent of the objects. Interestingly, they find that 
the fraction of young galaxies increases with decreasing 
galaxy mass and decreasing environmental density. They 
also find that the environmental dependence is most pro- 
nounced at low a and disappears at high a (see Fig. 8 
in their paper). Despite the considerable differences in 
sample selection and analysis, ou r results are r e mark - 
ably consistent with one another. iThomas et al.l (|2010f ) 
further show that the young galaxies have lower a/Fe ra- 
tios than average and most of them show signs of ongoing 
star formation through their emission line spectra, in line 



with our conclusion that recent star formation is likely 
the cause of the environmental dependence at low a. 

We emphasize that we have assumed that age is the 
dominant factor that determines the shape of the optical 
SED. The dependence on velocity dispersion and envi- 
ronment could also be caused by differences in metal- 
licity, differences in the contribution from blue horizon- 
tal branch stars or blue stragglers, or perhaps some- 
thing else. To accurately interpret the spectra is an 
extremely difficult task due t o the lack of calibration 
data (see the recent reviews bv lConrov eran[2boa .20101 
IConrov fc Gunnir2010l and references therein). 

Nevertheless, recent detailed studies show that star 
formation inde e d occurs in man y early- type g a laxies 
(Yi ct al.' '2005'; Rich ct al.' '2005'; 'Kavirai e t all [20071: 
j Schawinski ct al. 2007; Tcmi ct al.. 2009; Young et al l 
I200I iSalim fc Richll2010l) . iSalim fc RichI (|2010D, using 
high resolution far-UV (FUV) imaging with the Solar 
Blind Channel of the Advanced Camera for Surveys on- 
board the HST, show that for most early-type galaxies 
with recent star formation, it takes the fo rm of wide or 
conce ntric U V rings. The SAURON t eam (jShapiro et al.l 
I2OIOI see also iKuntschner et al.ll2010t ). find that star for- 
mation in early-type galaxies happens exclusively in fast- 
rotating systems and occurs in two distinct modes. In 
one mode, star formation is a diffuse process, correspond- 
ing to widespread young stellar populations and high 
molecular gas c ontent, possibly du e to (mostly minor) 
mergers (see also lKavirai et al.ll2009t ). In the other mode, 
star formation is concentrated into well-defined disc or 
ring morphologies, which may have been caused by re- 
juvenat ions within previously qu iescent stellar systems. 
Indeed, IKuntschner et al.l ()2010[ ) find that the most ex- 
treme cases of post-starburst early-type galaxies, with 
SSP-equivalent ages of < 3 Gyr, frequently show signs of 
residual star formation and are generally low-mass sys- 
tems. Therefore, variations in recent star formation is 
a plausible cause for both strong velocity dispersion de- 
pendence and weak but significant environmental depen- 
dence we see here. 

4.2. Lick indices 

The age-metallicity degeneracy has haunted stel- 
lar population analysis for decades. Nevertheless, a 
promising approach to breaking it remains the com- 
bined use of m ultiple absorption- line indices (e.g. , 
Faber et al.l [198 5: Gonzalez 1993; Worthev et all 11994 
Worthev fc Ottaviani.1997; Trager et all 119981 ). In this 
and the next section we measure the absorption-line in- 
dices of our spectra and compare them with state-of-the- 
art SSP models in order to infer their ages, metallicities, 
and a— enhancements. 

4.2.1. Lick indices of smoothed flux-calibrated spectra 

Lick indices measure the absorption-line strength of 
features in the SEDs and are widely used in studying stel- 
lar populations. The standard Lick indices are measured 
on low-resolution spectra that are observed with the Lick 
IDS instrument and not flux-calibrated. The SDSS did 
not observe the bright stars in the Lick library; therefore 
we cannot match the fluxing and calibration of the Lick 
system. Therefore we attempt an intermediate step by 
measuring the Lick indices of smoothed flux-calibrated 
SDSS spectra. 
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The bandpasses of t he Lick indices s tudied here are 
defined in Table 1 of IWorthev et alJ (|1994l ) and Ta- 
ble 1 of IWorthev fc Ottavianil (I1997D . To calculate the 
Lick indices, we first smooth each spectrum to a res- 
olution that is equivalent to a velocity dispersion of 
(T = 325 km s~"^, as observed under the SDSS resolution 
(69 km s~^), to match the largest a of our sample. We 
then calculate the Lick indices with a modified version 
of the Lick_EW routine in Ez _Ases'^ package developed 
bv lGraves fc SchiavonI (|2008D . The Lick_EW routine re- 
ports the err ors of each Lick inde x calculated in the way 
suggested bv lCardiel et all ()1998[ ). 

The combined resolution of ct = 325 km s^^ and the in- 
strumental resolution of the SDSS (69 km s^^) is at lower 
resolution than the Lick IDS instrument resolution at 
which the Lick indices are defined. We therefore need to 
correct the measurements back to the effective Lick IDS 
resolution. To do so, we use the correction factors given 
by S07 in their Table 46, which are for a 14.1 Gyr old stel- 
lar population with solar metallicity. The factors vary at 
the few percent level for stellar populations with differ- 
ent ages and metallicities; we ignore those tiny differences 
here. We also note that these empirical corrections fo r 
a may have large uncertainties (e.g., Tr ager et al.lll998l ). 
which unfortunately are hard to quantify and we also 
ignore here. 

In the end, the Lick indices measured in this way are 
matched in resolution, but still not with the fluxing of 
the Lick IDS system. S07 gives empirical zero-point cor- 
rections (Table 1 in their paper). When necessary, we 
adopt these corrections to compare flux-calibrated Lick 
indices with standard ones. These corrections have very 
large uncertainties; therefore we remind our readers that 
any comparison of our measurements with the standard 
Lick indices should be taken with some caution. 

4.2.2. Emission line infill correction 

One of the most challenging issues in Lick index mea- 
surements is emission line contamination. The Balmer 
absorption features are contaminated by emission from 
ionized gas, either from star formation, AGN activity, 
or interstellar shocks. This emission line infill severely 
affects the measurement of the Balmer indices, which 
are thought to be the best age indicators. For example, 
a 0.1 A contamination to the absorption- line equivalent 
width (EW) of Hz? from emission translates to an SSP 
age 1 Gyr older (e.g., S07). 

We correct for emission-line contamination using 
[O III] A5007, which is correlated with the Balmer emis- 
sion. Similar to the template fitting in §4.1.4, for each 
individual spectrum, we mask out the emission line re- 
gions and fit the continuum to a nonnegative linear com- 
bination of a set of SSP templates. Unlike in §4.1.4, 
where we only use models of old and young stellar pop- 
ulations, we here use more templates, including some of 
intermediate age as well. We then subtract the best-fit 
model and measure the flux and equivalent width (EW) 
of [O III] A5007, which is more reliable than the di- 
rect H/3 emission measurement. We then correct for the 
H/3 emissi on hne assuming AH/3 = 0.6 EW(fO III] A5007) 
(see, e.g. jGonzalelliggsTlTVager et al.ll2000bO . We also 



correct II7 and 11(5 for emission assuming the decrements 
H7/H/3 = 0.46 and HJ/H/? = 0.26. For reference, the 
median II/3 correction for all galaxies in the Elliptical 
sample is ~ 0.37 A, ~ 4 Gyr in SSP analysis. 

We emphasize that the corrections we adopt are 
by no means perfect. The relation between H/3 and 
[O III] A5007 has a large scatte r (see also, 
iTrager et al.' '2000b'; 'Mehle rt et all l2000t iNelan et al.l 
2005; Kuntschncr ct al. 20061). For EW([0 III] A5007)= 
0.5 A, an error of 0.2 in the conversion factor translates 
to 0.1 A in the H/3 measurement. On the other hand, 
the [O III] A5007 has been widely used in the literature 
and is reliable to measure. 

Another way to correct for the emission is to use 
Ha emission measurements, by adopting the intrinsic 
Balmer decreme nt Ha/H/3= 2 . 86 an d a reasonable dust 
extinction (e.g., iGraves et ahl I2007D . If dust is well- 
behaved, this relationship should be tighter than that 
for [O III] A5007. However, the Ha and H/3 emission 
lines are entangled with the underlying stellar absorp- 
tion feature. In fact, they are almost impossible to mea- 
sure reliably if they are relatively weak (< 1 A), which is 
the case for most of our spectra. In particular, because 
measuring the emission lines requires modeling the stel- 
lar continuum, using corrections based on Ha tends to 
merely recover the H/3 index of the continuum template 
itself. 

Here we estimate the emission line infill using 
[O III] A5007 because the line is easier to measure in- 
dependently of the stellar continuum, and because its 
intrinsic scatter is less important due to the large sam- 
ple size in our case. In any case, we have tried using 
Ha by adopting E(Ha-H/3)= 0.1 mag and find that the 
dependence of H/3 absorption-line strength on a and en- 
vironment is basically not affected. 

4.2.3. Velocity dispersion dependence of Lick indices 

Figures [13] and [14] present our Lick index measure- 
ments as a function of velocity dispersion and environ- 
ment. Figure [T3I shows the Balmer absorption lines: H/3, 
H7A, H7F, H(5a, and H(5f. Figure [Ml shows the metallic- 
ity indicators: (Fe), the average of Fe5270 and Fe5335, 
Mg5, [MgFe]', C24668, and Ca4227. The index [MgFe]' is 
defined as follows: 



* |http: //www.ucollck. org/'^graves/EZ_Ages .html| 



[MgFe]' = ^/Mgh- (0.72 • Fe5270 + 0.28 • Fe5335). 

[MgFe]' is a good metallicity indicator almost indepen- 
dent of a/Fe ratio variations (TMB). The left columns 
of these figures show the field galaxies and the middle 
columns show the rich-group galaxies. We single out the 
Coma galaxies in the middle column and show them sep- 
arately. We exclude the poor-group galaxies in these two 
columns, though we do show linear fits to their distribu- 
tion. We also show the indices measured on the average 
spectra in the right column. 

The most obvious feature from these plots is the strong 
index-cr relation, whi ch has been studied hitensively in 
numerous works fe.g..lTer levich et al.lll98ll: 'Gorgas et_a]J 
19901: iGuzman et al.lll99 2: Bender ct al. 1993; J0rgcnse3 
1997t [Be nder et al.' 1998; Bcrnardi ct al. 199? 
Colless et a l. 1999; J0rgcnsen 1999; Co ncann on et al] 
20001 [Ku ntschner 2000; Poggianti et al.l l200l[ 
Proctor fc Sansoin ,2002: Bernardi et al.l l2003b[ 
Caldwell et all lloM IWorthev fc Collobert 120031 



10 



Zhu et al. 



iMehlert et aH [200l iNelan et alJ [20051 T05). We per- 
form a linear least-squares fit to each distribution with 
the following form: 

Index = ci+C2 (logio ^ 2.20), (4) 

and list the results in Table [21 We show these fits in each 
panel of Figures [T3| and [T^j for all three group classifica- 
tions. 

All the Balmer indices strongly anti-correlate with a 
(c2 < 0). Their strength is weaker at higher ct, indicat- 
ing older ages or weaker recent star formation for more 
massive elliptical galaxies. Of these indices, H/3 is the 
one most investigators rely on as the best age indicator, 
despite the fact that it is the one with the most signif- 
icant emission line contamination issues. In particular, 
it is less s ensitive to metal lines than are other Balmer 
lines (e.g.,[KomeEaD[2005). 

The right column of Figure [T3| shows the measurements 
of the average spectra. All such measurements appear to 
follow the best-fit scaling relations very well except for 
the II/3 measurement in the lowest a bin. This discrep- 
ancy is possibly due to strong Balmer emission compared 
to the [O III] A5007, which makes the emission correction 
underestimated. 

In Figure [151 we take a closer look at the aver- 
age spectra within the wavelength range containing H/3, 
[O III] A5007 and Ha, where we also overplot the best- 
fit template to the continuum for the rich-group galax- 
ies. The Balmer emission lines are indeed very strong 
compared to [O III] A5007 in the average spectra of the 
lowest a bin, especially for the field galaxies. Specif- 
ically, direct measurements for the field galaxies yield 
EW(Ha)= 2.77 A and EW([0 III] A5007)=0.87 A. Us- 
ing [O III] A5007 gives the H/3 correction 0.52 A, while 
using Ha gives 0.88 A, larger by 0.36 A. 

These average spectra are each stacked from spectra 
of ~ 100 galaxies: thus, one galaxy with very peculiar 
emission lines can introduce such a feature in the aver- 
age spectra. In contrast, the fit to the scaling relation is 
less affected by a single data point outlier, so we will rely 
below extensively on the scaling relation. These results 
emphasize that achieving the emission line infill correc- 
tion at the required level of precision is extremely diffi- 
cult! 

Ah the metalhcity indicators, (Fe), Mg&, [MgFe]', 
C24668, and Ca4227, correlate with a (c2 > 0), though 
to varying degrees. The increasing strength with increas- 
ing a indicates higher metallicity for more massive ellip- 
tical galaxies. The tight relation with a of Mg6 (slope 
C2 ~ 3.44) is much stronger than that of (Fe) (c2 ^ 0.93), 
implying a higher [Mg/Fe] ratio for more massive galax- 
ies. Another a-element indicator C24668 also appears 
to be more strongly correlated with a than (Fe) with a 
slope C2 ~ 5.78. Taken together, these results are evi- 
dence of stronger a— enhancement in more massive ellip- 
tical galaxies. In fact, T05 show that increasing metal- 
licity, a— enhancement and older age account for the 
Mg- cr relation by 60%, 23% and 17% respectively (see 
also IMehlert et al.ll2003| ). The dependence of Ca4227 on 
a appears to be much weaker than other a-elements, 
with a slope C2 ~ 0.58, which implies that more raaB- 
sive galaxies may be mor e calcium-underabundant (see, 
e.g. iThomas et alll2003al) . 

To summarize, we see signs that the stellar popula- 



tions in more massive elliptical galaxies are older or with 
weaker recent star formation, more metal-rich, and more 
strongly a— enhanced, as has been found in the past. We 
will come back to this in the next sections when we in- 
terpret our index measurements using SSP models. 

4.2.4. Environmental dependence of Lick indices 

Figure [T31 shows that the Balmer index-cr relation for 
elliptical galaxies varies as a function of environment. 
The Balmer lines are systematically stronger (~ 0.1 A for 
H/3) in field galaxies than in rich-group galaxies. The 
difference vanishes at high a. This agrees well with the 
environmental dependence we see in the continua of the 
average spectra (Figures [TUl and ITU) . We also find that 
the Coma galaxies in our sample all have weaker Balmer 
lines than average. 

Balmer indices are thought to be the best age indica- 
tors because massive stars (e.g., A stars) in a young stel- 
lar population exhibit strong Balmer-line absorption in 
their integrated spectra. The environmental dependence 
thus suggests that the stellar populations in field galaxies 
are in general younger than those in group galaxies. 

In contrast, the distributions of metallicity indicators 
are remarkably similar for elliptical galaxies in different 
environments. The (Fe) index is nearly identical for all 
the samples. The a-element indices seem to be slightly 
stronger in rich-group galaxies, but only at a barely de- 
tecta ble level (as seen previously by, e.g.. lBernardi et all 
119981) . However, the same index measurement does not 
necessarily mean the same metallicity. At a given cr, if 
group galaxies are older than field galaxies as the Balmer 
indices indicate, then having the same metallicity index 
strength means they are more metal-poor. 

In this section, we studied the Lick indices, their 
relation with cr, and their environmental dependence. 
Consistent with previous studies, we have found strong 
index-cr relations, indicating higher metallicity, stronger 
a— enhancement and older ages for more massive ellipti- 
cal galaxies. We have also found that Balmer indices are 
systematically weaker in group galaxies, implying older 
ages for group galaxies. We do not see an obvious de- 
pendence on environment for (Fe). The a-element in- 
dices (Mg h, C24668 and Ca4227) show shghtly stronger 
absorption in group galaxies, but only at a barely de- 
tectable level. This result suggests that group galaxies 
are slightly more iron-poor and slightly more strongly 
a— enhanced than field galaxies. 

However, to convert the measurements to more quan- 
titative parameters of the stellar populations introduces 
the uncertainties associated with the stellar population 
synthesis. Later in the paper, we will nonetheless make 
an attempt to use state-of-the-art SSP models to study 
the stellar populations in elliptical galaxies. 

4.2.5. Systematic Uncertainties 

We conclude this subsection by summarizing the sys- 
tematic uncertainties affecting our Lick index measure- 
ments. 

We measure the Lick indices on the flux-calibrated 
spectra that are smoothed to cr = 325 km s~^ under 
the SDSS resolution 69 km s~^. We have used empiri- 
cal corrections to correct the indices back to the effec- 
tive Lick IDS resolution. However, the corrections may 
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have large uncertainties (jTrager et al.lll998[) . The stan- 
dard Lick indices are measured on spectra that are not 
flux-calibrated. Therefore, any direct comparison of our 
measurements with standard ones should be taken with 
caution. 

The SDSS spectra are taken with fibers that enclose 
an aperture with a diameter of 3". Since we limited our 
sample to be lower than redshift 0.05, the spectra are of 
the central region of the galaxies. The velocity disper- 
sion is therefore essentially the velocity dispersion in the 
central region. For galaxies of different angular sizes, the 
aperture therefore encloses the central regions of differ- 
ent sizes. Recent work has shown that the radial profile 
of a follows a power law or oc W w here R is radius 
and the index 7 ^ ~P-04 to —0.07 (e.g., |j 0rgense n et al.l 
[19951 IMehlert et al.l[200l ICappellarietal. .20061 The 
effective radii of the Elliptical sample have a median 
value ~ 7". If we correct all the cr to that at effective half 
light radius with such a scaling relation, this translates 
to a difference of ^ 3%; this small effect should not affect 
the analysis presented here. 

The Lick indices measured here also come from the 
spectra of stellar populations in the central regions, with 
a much smaller fraction from those in the outer region 
within the aperture. The Lick indices also probably cor- 
relate with radius, and the gradients of indices such as 
Mg h and (Fe) seem to correlate with that of a as well. 
On the other hand, H/3 is roughly constan t with ra- 
dius (e.g., Mchlcrt ct al. 2003; Kuntschncr e£aD 120061 : 
ISanchez-Blazquez et all i2006ait iRawle et al.i i2010l ) . The 
correlations are still poorly understood and we do not 
correct for them. 

Finally, the emission line infill correction for the 
Balmer indices of an individual spectrum is extremely 
difficult, even in the case of stacked spectra with ~ 100 
galaxies — one galaxy with a peculiar emission line ra- 
tio can introduce strong errors in the correction. Nev- 
ertheless, we are using [O III] A5007 to correct Balmer 
emission lines, which should be a good compromise for 
such a large sample. We also have tried using Ha too, 
which basically does not affect the relative dependence 
on cr and environment. 

All of these factors may introduce some errors at the 
few percent level in our measurements. We do not expect 
they introduce any significant systematic bias in the dif- 
ferential aspect of our analysis, because they should not 
correlate strongly with which subsample we consider. 

4.3. Comparison of Lick indices with SSP models 
4.3.1. SSP models 

One purpose of stellar population synthesis models is 
to place constraints on the history of star formation and 
chemical enrichment of galaxies from their integrated 
SEDs. We apply two state-of-the-art versions of these 
models to our measurements: those of TMB and S07. 

There are a few caveats that make our goal very chal- 
lenging. For example, age and metallicity ha ve extremely 
simil a r effects on the c olor and SEDs (e.g.. lFabeilll972L 



1973; 



198a 



Oconnelll ll98C)t iRosd 119851 : iRenzini fc Buzzonil 
Worthev et a'iTll994D : Lick indices have very 



l2005h : finally, many elliptical galaxies have non-solar 
[g/Fe] w hich further compli cates t he modeling (e.g. , 
iPeletiei^ n j8l IWorthev et all 119921: iDavies et all 119931: 
;McWilliam"XRic31l99J^ 

Recent SSP models have incorporated adjustable 
abundance patterns for multiple elements and allow for 
more reliable derivation of age, metallicity, and element 
ratio from absorpti on-line indices (e.g., 'Borges et al 



19951 IWeiss et al.lll99 5: Tantalo ct al. 1998; Tragcr.,eial 
2000al. TMB. iMarastonI 120051 : iCoelho et al.ll2007L S07) 



broadly-defined line windows which makes dire ct trans 
lation into element abundances difficult (e.g.. iGreggic 
Il997t iTantalo et anil998l: iKorn et al.|[2005t IServen et al 



Most of these models, however, are based on the Lick 
system and are not directly applicable to our data, as dis- 
cussed in §4.2.11 Nonetheless, S07 provide empirical zero- 
point corrections to correct flux-calibrated Lick indices 
onto the standard Lick system (see their Table 1). We 
apply these corrections to our Lick index measurements 
when we compare them with the TMB models. The S07 
models, meanwhile, are built up on flux-calibrated spec- 
tral library, so no corrections are necessary for that com- 
parison. 

The corrections necessary for the comparison to TMB 
have a very large scatter (see Table 1 and Figure 1 of 
S07). Therefore, while we will compare to TMB below, 
our main focus will be on the models of S07. 

4.3.2. Comparison with TMB models and S07 models 

Before deriving the stellar population parameters of 
our sample with SSP models, we first compare our mea- 
surements to the predictions of each model on grids of 
age, metallicity and a— enhancement. 

One of the main differences between TMB models and 
S07 models is that TMB models are built at fixed metal- 
licity [Z/H], while S07 models are at fixed iron abundance 
[Fe/H]. To convert between [Z/H] and [Fe/H], we adopt 
the relation given by TMB: fz/H] = [F/H]-h0.94 [a/Fe]. 
Because we tie other a-elements to [Mg/Fe] when build- 
ing up S07 models (see next subsection), we assume 
[Mg/Fe] = [Q;/Fe] in the comparison. 

In Figure [16] and [iTl we compare our measurements 
of H/3 versus (Fe) with the S07 and TMB models as- 
suming [a/Fe] = [Mg/Fe]=0.3. We choose [a/Fe]=0.3 be- 
cause most elliptical galaxies have super-solar [a/Fe]. 
In the following analysis we only use elliptical galaxies 
with cr between 125 km s^^ and 200 km s"^ (2.10 < 
logj^Q (T < 2.35, the middle a bin), to ensure that we 
do not confuse cr dependence with environmental de- 
pendence. The median a values within this range are 
166 km s~^, 167 km s^^ and 172 km for field, poor- 
group, and rich-group galaxies, respectively. We also 
show the median measurements of each subsample as 
filled stars and the la scatters as error bars below the 
legend. Finally, for comparison, in the lower right panel 
of each figure we compare the physical parameters of the 
galaxies in our sample to those of the early-type galaxy 
sample of T05. 

In general, the S07 an d T MB models agree well ( see 
S07; [Graves fc Schiavonll2008l : iKuntschner et al.l l2010 for 
more thorough comparisons between the two models). 
We see that the galaxies in our Elliptical sample have 
a median SSP-equivalent age of 7 Gyr. These measure- 
ments show that group galaxies, especially rich-group 
galaxies, appear to occupy a different locus than field 
galaxies. In particular, the distribution of group galaxies 
extends to older ages and lower iron abundance: rela- 
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tive to the other samples, the rich-group galaxies ap- 
pear to more commonly have sub-solar iron abundance 
([Fe/H]< 0) and older age (> 7.0 Gyr); meanwhile, the 
field galaxies are typically < 7.0 Gyr. This result sug- 
gests that the rich-group galaxies are slightly older and 
slightly more iron-poor (in terms of [Fe/H]). We also 
show the T05 sample in the bottom right panel and see 
very similar effects. Note when comparing with S07 mod- 
els, we have corrected measurements of the T05 sample 
to the flux-calibrated indices. We also see that the Coma 
galaxies are systematically older than average. 

The strength of the (Fe) index is mainly determined 
by the iron abundance. For total metallicity, [MgFe]' is 
a better indicator. In Figure [T51 and [T51 we plot our mea- 
surements H/3 vs. [MgFe]' against S07 and TMB mod- 
els with [a/Fe]=0.3. These measurements show that in 
different environments the total metallicities [Z/H] are 
very similar. The T05 sample^, however, exhibit slightly 
lower total metallicity for ellipticals in high-density en- 
vironments. 

In Figures [201 and we compare our measurements of 
(Fe) and Mg b against S07 and TMB models with an age 
of 7 Gyr. We see that almost all elliptical galaxies have 
super-solar [Mg/Fe], with a median value close to 0.3, 
and the elliptical galaxies in rich groups exhibiting some- 
what higher [Mg/Fe]. In the lower right panel, we also 
show the elliptical galaxies in the T05 sample. Although 
T05 do not claim to see an environmental dependence 
in their analysis, which includes lenticular galaxies (not 
shown here) , we see the elliptical galaxies in high-density 
environments in their sample also appear to be slightly 
more strongly a— enhanced. 

4.3.3. SSP- equivalent parameters: age, [Fe/H], [Mg/Fe] and 

[Z/H] 

Here, we use S07 models to derive the SSP-equivalent 
parameters. To do so, we create a grid of models and 
interpolate our results onto the grid. 

First we use the deltabund code in the EZ_ages pack- 
age to create SSP models on a grid of three parame- 
ters: age, [Fe/H] and [Mg/Fe]. We use the solar-scaled 
isochro nes as suggested in the EZ_ages documentation, 
and the 'Salpeterl (fl955') IMF. We set [0/Fe] to be zero 
for these isochrones and tie other a-elements to Mg. The 
age allowed in the models ranges from 1.2 Gyr to 17.7 
Gyr, the [Fe/H] ranges from —1.3 to 0.2 and we generate 
models at [Mg/Fe] between —0.1 and 0.5. W e refer the 
readers to S07 and [Graves fc Sch iavon' (2008) for details 
about S07 models and the EZ_ages package. 

Once we have all the models, we employ a two-step 
interpolation method to derive age, [Fe/H] and [Mg/Fe]. 
For each galaxy, we first calculate age and [Fe/H] using 
H/3 and (Fe), at all [Mg/Fe] (from -0.1 to 0.5). We then 
compare the pair (Fe) and Mg h with models at the me- 
dian age we obtained in the first step to calculate a new 
[Fe/H] and [Mg/Fe]. We then update the age by inter- 
polating the ages found in the first step at the [Mg/Fe] 
we found in the second step and iterate the second step 
with the new age. The iteration usually only needs ~ 2 
steps before fully converging. After convergence, we have 

^ The public catalog of T05 only includes the average index 
(Fe), so we show instead the combination [MgFe] = -^Mg b ■ (Fe). 
However, the difference is very small. 



the parameters age, [Fe/H] and [Mg/Fe] for each galaxy 
with H/?, (Fe) and Mg6. We also give [Z/H] assuming 
[Z/H] = [F/H]-K0.94[a/Fe]. For those measurements that 
fall out of the model grids, we set the parameter to be the 
boundaries, i.e., the maximum or the minimum in the 
models. We have also tried using higher-order Balmer 
lines instead of H/3 and found very consistent results. 

Figure [22] shows all the measurements as a function of 
environment and a. We also overplot the derived param- 
eters for the average spectra as large open symbols. In 
addition, we also show as large filled symbols the derived 
parameters at the same a of the average spectra but as- 
suming that the index-cr relations (Eq. 1 and Table 2) 
hold exactly. 

As we pointed out above, the H/3 absorption in the av- 
erage spectra is underestimated due to the underestima- 
tion of the emission correction. The derived age therefore 
is overestimated, and the derived [Fe/H] underestimated. 
For this reason, in this context we trust results based on 
the index-cr relationships more than the average spectra. 

The most obvious feature in Figure [H] is the strong 
dependence on a of all three parameters, as we expect 
from the index-cr relations (though with a large scatter 
due to the combined errors from the three indices). More 
massive elliptical galaxies are older, more metal-rich and 
more strongly a— enhanced. The environmental effect 
on age is also apparent: galaxies in groups appear to be 
older than their counterparts in the field. Meanwhile, 
we do not see an obvious environmental effect on [Fe/H], 
[Mg/Fe] and [Z/H], implying the environmental depen- 
dence of [Fc/H], [Mg/Fe] and [Z/H] is very subtle, if it 
exists. 

To quantify the a and environmental dependence of 
the SSP-equivalent parameters we fit the scaling relation 
for each distribution as a function of cr as follows: 

SSP parameters — si+ S2 (log^o f — 2.20). (5) 

However, we need error measurements for the derived 
parameters in the fitting. The errors are unfortunately 
more difficult to quantify. 

We assume all the measurement errors are independent 
and propagate the errors to the derived parameters with 
a Monte Carlo method. To do this, we create 50 fake 
measurements for each galaxy by adding Gaussian errors 
to the measurements of H/3, (Fe) and Mg &, and derive 
the new parameters for the fake measurements. We then 
calculate the standard deviations of the parameters of 
the Monte Carlo ensemble for each galaxy and quote it 
as the Monte Carlo error of the derived parameters. 

We show the median Monte Carlo error of each pa- 
rameter in the upper left corner in Figure [22l However, 
the individual error of each point can be heavily biased, 
due to edge effects on the grid. Almost all the galax- 
ies have some simulated data points falling outside the 
grid. This effect is more severe for those galaxies with 
parameters closer to the boundaries. In the extreme case, 
the simulated measurements of a galaxy with index mea- 
surements lying outside the grid are almost always also 
outside the grid; thus, for the Monte Carlo ensemble of 
such a galaxy, the parameters derived with our method 
are almost always the same as the boundary values in 
the models. Therefore, the Monte Carlo errors are unre- 
alistically small for the parameters of the galaxies with 
measurements close to or outside the model boundaries. 
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However, for most of the derived parameters that are 
relatively far away from the boundaries, the Monte Carlo 
errors are very similar. We therefore apply a uniform er- 
ror (the median Monte Carlo error) to each parameter 
of all galaxies. We also ignore the galaxies with mea- 
surements that are out of grids when fitting the scaling 
relation, because they have the same values for the de- 
rived parameters (either the maximum or the minimum) 
at all (7, which will heavily skew the fitting. 

Using these choices, we fit equation ([S]) to the derived 
parameters. Figure[22l shows the best-fit scaling relations 
and Table [3] lists all of the coefficients. Besides the errors 
in the slopes and intercepts, we have also calculated the 
internal scatter of the distribution of the whole sample 
in our three cr bins. The internal scatter (in dex) for 
log;^o s-gs, [Fe/H] and [Mg/Fe] in the three a bins (in as- 
cending order) are: (0.22, 0.21, 0.18), (0.24, 0.14, 0.10) 
and (0.13, 0.11, 0.08), respectively. Recall that the me- 
dian a in the three bins are 100, 167 and 232 km s~ , 
respectively. 

The strong a dependence is obvious from the best-fit 
relations. Elliptical galaxies at the highest masses are 
older, more metal-rich, and more strongly a— enhanced. 
The rich-group galaxies are also apparently older than 
the field galaxies, by ~ 1 Gyr. We note, however, that 
the difference between different subsamples is relatively 
small compared with the intrinsic scatter. In the previous 
sections, we showed that the environmental dependence 
of the average spectra and H/3-cr relations vanishes at 
higher a. We therefore expect that the SSP-equivalent 
age dependence on environment may correlate with a 
too. Indeed, the slope of the age-a relation for field galax- 
ies is slightly larger than that for rich-group galaxies, but 
the difference is very small. The dependence of the age- 
environment relationship on a thus may be buried in the 
large scatter in this diagram. 

Aside from the strong a dependence and the obvious 
environmental dependence of age in the best-fit scaling 
relations, we do see very subtle environmental effects on 
[Fe/H] and [Mg/Fe]. However, these are at a barely sig- 
nificant level (Icr). And we do not see obvious depen- 
dence of total metallicity [Z/H] on environment. Com- 
pared to the field galaxies, the rich-group galaxies are 
slightly more iron-poor (only by ^ 0.01 dex in terms 
of [Fe/H]) and slightly more strongly a— enhanced (only 
by ~ 0.01 dex). These differences are consistent with 
the environmental dependence of the index-tr relations 
we saw in §4.2.41 In Figure [T4] and Table [21 we show 
that the (Fe) absorption-line strength is almost identical 
in different environments. However, because field galax- 
ies are younger, they need slightly more iron to produce 
the same (Fe) signature as the rich-group galaxies; field 
galaxies also have slightly weaker Mg 6 (and other a- 
element indices), implying slightly lower a— abundance 
in the field, as we see here. 

The environmental dependence we see agrees very well 
with previous studies (e.g., T05; iBernardi etal] 120061: 
ICooper et aLll2010t ). T05 show that early- type (E/SO) 
galaxies in high-density environments (mostly in the 
Virgo and the Coma clusters) are older by ~ 2 Gyr 
than galaxies in low-density environments. We have also 
shown that the Coma galaxies in our Elliptical sam- 
ple are systematically older than average; thus, the Coma 



galaxies in their sample may have contributed to the rel- 
atively stronger dependence in their results. They also 
found a relatively stronger dependence for metallicity; in 
particular, they found that galaxies in high-density en- 
vironments are more metal-poor by 0.05 — 0.1 dex (in 
terms of [Z/H]) with no dependence for [a/Fe]. As seen 
in Figure [201 if we only look at the elliptical galaxies in 
their sample, we see the galaxies in high-density environ- 
ments appear to be slightly more strongly a— enhanced. 
Considering the large scatter in the SSP parameter-cr re- 
lation, our results agree very well with each other. 

By compari n g sa mples in different environments, 
IBernardi et al.l ()2006l ) find that the spectroscopic differ- 
ences between early-type galaxies in high-density envi- 
ronments and in low-density environments are very sim- 
ilar to the differences between early-types at redshift 
z — 0.17 and z = 0.06. Under the reasonable assumption 
that the primary difference between the samples at each 
redshift is overall age, they therefore conclude that early- 
type galaxies in high-density environments are 1 Gyr 
older than those in low-density regions. Thus, our con- 
clusions are in very good agreement with theirs. They 
also do not find noticeable environmental dependence in 
total meta llicity. Combining res ults from different groups 
(e.g., T05. IBernardi et al.|[200l and this work), we note 
that the environmental dependence in metallicity is in- 
deed very subtle. 

Usin g samples drawn from the SDSS, ICooper et al.l 
remove the mean dependence of average overden- 
sity (i.e., environment) on color and luminosity, and find 
that there remains a strong residual trend between stellar 
age and environment, such that galaxies with older stellar 
populations favor regions of higher overdensity relative to 
galaxies of like color and luminosity. Their conclusions 
therefore are also consistent with the environmental de- 
p endence we find here . 

iThomas et al.l (|2010( l do not find a noticeable environ- 
mental dependence for the bulk old population of early- 
type galaxies, likely because they separate the young ob- 
jects with signs of recent star formation from the bulk 
old population. In §4.1.4! we argued that the environ- 
mental dependence of the average spectra was likely due 
to recent low-level star formation in low a galaxies. We 
conc lude, therefore , that our results would be consistent 
with [Thomas et al] ([2010( ) if they had included the young 
early-type galaxies in their scaling-relation analysis. 

In Figure 1231 we compare the linear relationships we 
derive between velocity dispersion and age, iron abun- 
dance, Qf-enhancement, a nd total metallici ty against 
those derived by T05 and iNelan et al] ([20051). We em- 
phasiz e that the fits derived by both T05 and [Nelan et al.l 
(2005) include elliptical and lenticular galaxies, while our 
sample only inc ludes elliptic a l gala xies. We determine 
the intercept of [Nelan et al.[ ([2005| ) from their Table 8 
and Figure 13. For T05, we use the mean of the slopes 
and the intercepts for samples in high-density and low- 
density environments. We also overplot the Icr scatter in 
the distribution of the parameters for the whole sample 
as error bars (though of course the errors in the mean 
values are much smaller). 

Figure [23[ shows that our slopes are in very good agree- 
ment with theirs. Considering the large scatter in the 
distribution, uncertainties in zero-point corrections, dif- 
ferent samples and measurements used, and different 
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models used, we find the agreement is remarkable and 
very encouraging for SSP analysis. The main difference 
is that our age-cr relation is systematically younger, by 
~ 3 Gyr. This is mainly caused by the difference in the 
H/3 measurements in that our H/3 indices are systemati- 
cally stronger than theirs, by 0.3 A (see, e.g.. Figure 

[ISD. 

The discrepancy in age can probably be attributed to 
the emission li ne infill corre c tion a nd / or sample bias. We 
notice that in iNelan et al.l (|2005[ ) they exclude galaxies 
with significant emission lines, by which they may have 
thrown away galaxies with significant young stellar pop- 
ulations. For the resulting sample without significant 
emission lines , they do not correct for emission-line con- 
tamination. In the sample compiled by T05, we find 
in the literature that the indic es of some galaxies are 
corr ected for emission ( |Gonzalcz l|1995 ), whi le others are 
not ()Beuing et al.ll2002llMehlert et al.ll2000t ^°). Note the 
median correction for H/? in our Elliptical sample is 
~ 0.37 A. Although this difference is fairly small com- 
pared to the absorption strength (~ 2 A), it translates 
to ^ 4 Gyr in SSP analysis. In addition, in T05 more 
than half of the galaxies in high-density environments 
are from the Coma cluster. As we have shown here, 
the Coma galaxies in our sample are systematically older 
than average. Taking into account the emission line in- 
fill correction, sample bias, and the large scatter of the 
distribution, we conclude that these age— cr relationships 
agree fairly well with each other. 

4.3.4. Caveats on SSP parameters 

We emphasize that the parameters calculated above 
are SSP-equivalent parameters. That is, they can only 
be taken literally if elliptical galaxies formed in a short 
enough, and uniform enough, burst of star formation. 
Such a scenario is probably incorrect. Recent obser- 
vations from GALEX, Spitzer, and HST have shown 
that a significant fraction of early-type galaxies exhibit 
strong UV excess, PAH emission and IR e xcess, implying 
possible low-leve l recent star formation (|Yi et al.ir200E 



Rich 



et all l2005t jKavirai et al.' '2007"; 'Schawinsk i et al 



2007 : lTemi et al.l [l)09; Youn g ct al. 2009; Salim fc Rich 
2010). Similarly, in HA.lAl we showed that the en- 



vironmental dependence of our frosting model param- 
eters (e.g., iTrager et all l2000al: iGebhardt et"aLl l2003l 
S07) w ere consistent with the results of lSchawinski et al.l 
(|2007( ) who used GALEX NUV photometry to look for 
recent star formation in early-type galaxies. 

If recent star formation is responsible for some or all of 
the (T and environmental dependence of the average spec- 
tra, then instead of giving the age of the dominant old 
stellar population, the derived SSP-equivalent age will be 
significantly younger (than the real age of the old popu- 
lation) because of t he existence of a sig nificant fraction 
of yo u ng stars (e.g. iTraeer et al.l[2000al : ISerra fc Traged 
I2007D . ISerra fc Trageil (|2007t ) show that the Balmer in- 
dices of a composite stellar population depends primarily 
on the mass fraction and age of the younger component. 
If so, then the a and environmental dependence of the de- 
rived age we see here may refiect partly, or totally, the a 
and environmental dependence of the mass fraction and 

We are not sure whether or not they correct emission for the 
galaxies with strong emission {EW(H/3)> 0.3 A). 



age of the young component in the frosting model. That 
is, elliptical galaxies at different tr may host the old base 
stellar populations of the same (or similar) age, but the 
SSP-equivalent ages are younger at lower a and in the 
field due to stronger recent star formation. 

Meanwhile, because of the age-metallicity degeneracy, 
the derived metallicity might be overestimated due to 
the underestimated age — assuming the metallicity in- 
dices are primarily tracing the old population. Precisely 
how the metal absorption indices behave in the case of 
composite stellar populations is as yet unclear. 

5. DISCUSSION 

In the above sections, we have found systematic differ- 
ences in the average spectra of ellipticals, as a function 
of a and environment. In addition, we have measured 
the dependence of SSP-equivalent ages, [Fe/H], [Mg/Fe] 
and [Z/H] on these parameters. 

As previous studies have shown, stellar ages of elliptical 
galaxy seem to be older at higher masses. In this sense, 
our results fit into the n ow-classic "downsizing" scenario 
that ICowie et al.l (|1996D presented in the context of the 
more general galaxy population. This age dependence 
on a appears to be explainable in the hierarchical sce- 
nario of galaxy formation, as l ong as star formation can 
be shut down effectively (e.g.. iDe Lucia et al.ll2006l ). In 
these models, the higher mass halos are comprised of 
small systems that on average collapsed earlier. If there 
is a mechanism to shut off their star formation quickly 
enough, these progenitors become gas poor and thus the 
higher mass systems will have older stellar population 
ages. 

The precise nature of what shuts off star formation in 
elliptical galaxies is unknown. AGN are a possible source 
for such feedback, t hrough b oth radiative and mech anical 
heating processes (iCiotti fc pstrikCT_1997; Benson et al.l 
20031 IDe Lucia et al.ll2006t IVernaleo fc Revnoldsl 120061 : 
Ciotti et al.l l2009f ). A commonly invoked mechanism 
for triggering the AGN is by feeding the central 
black holes during a major merger (e.g., 'S ilk fc ReesI 
1998: Mattco et al. 2005; Hopkins ct al. 2008). Indeed, 



Hopkins et all f2008i) argue that in these models the time 



since the last major gas-rich merger is an increasing func- 
tion of galaxy mass, and could explain an age-cr relation- 
ship of ellipticals. Other mechanisms for shutting off star 
formation also are potentially important: Ciotti ct al] 
((200l fuel AGN feedback using the rec ycled gas from 
dying stars, while iJohansson et al.l (|2009l ) and others in- 
voke gravitational heating from infalling stellar systems. 

Elliptical galaxies also have a strong dependence of 
a— abundance on a. This trend implies that the na- 
ture of chemical enrichment depends on galaxy mass. 
The iron-peak elements mainly come from Type la su- 
pernovae, while the a-elements are mainly produced in 
massive stars experiencing Type II supernovae. That 
the a— enhancement (or probably more accurately, the 
iron-deficit) is stronger in more massive elliptical galax- 
ies may imply that their star forming timescale is shorter 
than less massive elliptical galaxies, before the de- 
layed Type la supernovae enrich the star forming re- 
gions with iron-peak elem ents fe.g. lThomas "erallll998l: 
jPipi no fc Matteuccill200l T05. But also see lSmith et"all 
120091). Hierarchical mod els with feedback, such as those 
of iDe Lucia et al"] (|2006t ). appear consistent with this re- 
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suit at least qualitatively. 

Finally, ellipticals have a dependence of metallicity on 
stellar mass. If most of the stars are formed in situ, 
this trend might be explained by their deeper potential 
well, boosting chemical recyclin g by hampering the out- 
flow of gas from the galaxy fe.g.. lArimoto fc YoshiilHoSTl : 
iMatteuccH 1X9941) . However, if most of the stellar mass 
in ellipticals is produced in progenitor spiral disks, this 
metallicity trend may instead be a remna nt of the mass- 
metallicity trends k i iown for spirals (e.g. iGarnettI 120021 : 
iPilvugin et alll2004l: iTremonti et al.ll2004[) . 

Where the stars formed tha t now constitute ellipti- 
cal galaxies remains in debate. iNaab fc Ostriked (|2009f ) 
claim that the metallicities of present day spirals are 
several times lower than those measured for ellipticals, 
indicating that giant ellipticals cannot have b e en pro - 
duced from mergers of spirals. iHopkins et all (|2009bf ) 
counter that this problem could be mitigated if the stel- 
lar metallicity gradients of ellipticals are large enough; 
typical metallicity measurements of ellipticals only cover 
apertures that are a fraction of the effective radius, and 
the central regions might have highly enhanced metal- 
licities. It remains to be determined what would be 
the mass-metallicity relation predicted by a fully cos- 
mological simulation of the merger scenario (though see 
[Koba vashi 2005; Finlator & Dave 2008). 

A final wrinkle in the dependence of elliptical proper- 
ties as a function of mass is that their strength and scat- 
ter put some constraint on their late-time merger history. 
After all, if the largest galaxies are built from progenitor 
populations of lower-cr galaxies, one expects the depen- 
dence of age and chemical abundance on a to weaken 
with time. If the merging process is important to galaxy 
growth, these trends should be steeper and tighter in the 
past. An interesting quantitative prediction from the 
theoretical models would be how large this evolution re- 
ally was expected to be, given various merging scenarios. 

Alternatively, a variation of the stellar initial 
mass function (IMF) with galaxy mass — within 
the ellipticals or their progenitors — can ac- 
count fo r the metallic i ty and a— abundance trends 
as well (ICeiiia.rro et all l2003t iNagashima et al.l 120051: 
iKoppen et al.ll2007l : iHoversten fc Glazebrookll2008l ). 

Our results demonstrate that whatever processes pro- 
duce these correlations differ across environment, though 
only slightly. Explaining the differences between ellip- 
ticals in rich groups and in the field requires either a 
slightly later time of formation for field galaxies, or equiv- 
alently and more likely a slightly thicker frosting of recent 
star formation. Furthermore, these differences are sub- 
stantially stronger for low-mass galaxies than for high- 
mass galaxies. This distinction between the effects on 
high- and low-mass ellipticals is apparent both in the 
stellar continuum shapes, and in the age measurements 
based on the Lick index scaling relations (the age mea- 
surements we trust the most; large filled symbols in Fig- 
ure [22|). 

As noted above, the stellar ages of ellipticals may be 
relate d to the f o rmati on times o f their progenitor ha- 
los. iGao et al l (|2005D (see also IWechsler etldl 120061 : 
IZhu et all 120061 ) show that CDM dark matter halos in 
dense environments were assembled earlier than aver- 
age. Such an e ffect may also be reflected in Fig. 1 of 
iDe Lucia et al.l (|2006l ) , of the stellar age as a function of 



environment. If the ages of the populations are related 
to the formation time of the progenitor halos, these theo- 
retical results could explain the younger populations seen 
in the field. The theoretical results also predict that this 
effect declines at higher masses (reversing, in fact, above 
the current nonlinear mass). These results could be re- 
lated to the environmental dependence of mean stellar 
age. 

Alternatively, field ellipticals may simply have some- 
what better access to reservoirs of cold inflowing gas 
whose virial temperature in groups would be too hot to 
accrete at late times. In this case, the mean stellar age 
difference would derive mainly from recent star forma- 
tion, which as we note abo ve is actually seen in nearby 
fast-rotator ellipticals (e.g.. iKavirai et aLll2007t ). How- 
ever, it is unclear whether and why such an effect should 
be weaker for more massive ellipticals, as required by 
the dat a. For example, the gravi tational heating sce- 
nario of lKhochfar fc Ostriked (|2008l ) appears to predict a 
fairly strong dependence of age on environment for mas- 
sive galaxies. These considerations possibly favor inter- 
nal processes, rather than external ones, for regulating 
the star formation rates in the most massive ellipticals. 

The metallicity and a-abundance of stellar popula- 
tions are a much weaker function of environment than 
age. Whatever determines the chemical evolution of such 
galaxies must therefore be a weak function of environ- 
ment. For example, if the mass-metallicity relation is 
determined by a systematic variation in the IMF, our 
results would imply that the IMF does not vary sub- 
stantially with environment. In contrast, if merger his- 
tory is important in establishing the present day mass- 
metallicity relation, then ellipticals in different environ- 
ments must have similar enough merger histories. On 
the other hand, if the environmental dependence of age 
is caused by recent star formation ( §4.1.4p . we expect 
that stronger recent star formation would cause lower 
[a/Fe] and higher metallicity in field ellipticals. For ex- 
ample, by separ ating the young early - types from the bulk 
old population, iThomas et al.l (|2010') find that the young 
population have lower [a/Fe] and higher metallicity. It 
is possible that our sample is still too small to detect 
a stronger environmental dependence. It is therefore of 
great interest to investigate the relation of metallicity 
and a-abundance with environment with a larger sample 
of elliptical galaxies. 

6. SUMMARY 

We reanalyzed the images of red-sequence galaxies in 
the local universe using the homogeneous data set of the 
SDSS. By carefully examining the surface brightness pro- 
files, we selected 1,923 elliptical galaxies with velocity 
dispersion cr > 70 km s^^ at redshift z < 0.05. We found 
that elliptical galaxies dominate the bright /massive end 
(> L*) of the red sequence in the color- magnitude di- 
agram, and disk-dominated galaxies dominate at lower 
luminosity (< L*). We also found that elliptical galaxies 
and disk-dominated galaxies form different loci in color 
gradient- magnitude/velocity dispersion space, suggest- 
ing color gradient can be used for morphological clas- 
sification. 

We studied the dependence of properties of elliptical 
galaxies on velocity dispersion and environment. Group 
galaxies tend to have higher velocity dispersion, and thus 
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higher mass. 

We have calculated the average optical spectra as a 
function of velocity dispersion and environment, which 
all show a typical SED of an old stellar population. We 
found that the average spectra depend strongly on veloc- 
ity dispersion. Elliptical galaxies at lower-cr have a bluer 
continuum, stronger Balmer and nebular emission, and 
weaker metal absorption. This result is consistent with 
the color- magnitude/mass relation that brighter/more 
massive elliptical galaxies are redder in optical broad- 
band colors. 

Interestingly, we found weak but significant environ- 
mental dependence of the average spectra. Elliptical 
galaxies in the field have a bluer continuum, especially 
at wavelengths < 4000 A, and have stronger (but still 
weak) emission lines than their counterparts in groups. 
However, this dependence on environment appears pri- 
marily for low-cr ellipticals; the highest-cr ellipticals are 
much less affected. Assuming that age is the dominant 
factor shaping SEDs and that elliptical galaxies consist 
of an old base stellar population and a small frosting of 
young stars, we fit the frosting model to each spectrum 
in the Elliptical sample. Assuming these trends to be 
due to a young population, we found that the fraction 
of galaxies with significant young population is higher at 
lower velocity dispersion and in the field. We also found 
that the environmental dependence appears primarily at 
low velocity dispersion and vanishes at high velocity dis- 
persion. 

We measured the Lick indices of the flux-calibrated 
SDSS spectra of galaxies in our Elliptical sample. In 
agreement with previous work, we found strong index- 
<T relations. The Balmer absorption indices are stronger 
at lower velocity dispersion, while the metallicity indices 
are stronger at higher velocity dispersion. The Balmer 
indices of field galaxies are systematically stronger than 
those of rich-group galaxies. As in the case of average 
spectra, this dependence is most pronounced at low ve- 
locity dispersion and disappears at high velocity disper- 
sion. We did not find significant environmental depen- 
dence of metallicity indicators, only that the a-element 
indices appear to be slightly stronger in rich-group galax- 
ies. We also noted that the emission line infill correction 
for Balmer lines is an extremely difficult task. 

Assuming SSP analysis applies to elliptical galaxies, we 
have derived the SSP-equivalent age, iron abundance (in 
terms of [Fc/H]), a— enhancement ([Mg/Fe]) and total 
metaUicity ([Z/H]). We found that the SSP-equivalent 
parameters strongly correlate with velocity dispersion. 
More massive elliptical galaxies are older, more metal- 
rich, and more strongly a— enhanced. We also found 
that, galaxies in rich groups are systematically older than 
their counterparts in the field, by 1 Gyr. However, 
this effect is strongest at low velocity dispersion. We 
found that galaxies in rich groups are slightly more iron- 
poor and slightly more strongly a— enhanced, but only 
at a barely significant level. And we do not find no- 
ticeable environmental dependence of total metallicity as 

|http: //bias . cosmo . f as .nyu. edu/galevolution/elllptlcal| 



well. We have performed fits to the parameter-cr relation 
and found that our fits are in very good agreement with 
previous work, especially taking into account many un- 
certainties such as zero-point corrections, different sam- 
ples and data, and different SSP models employed. We 
found that emission line infill corrections can affect the 
age determination by ~ 3 — 4 Gyr. We also caution that 
the SSP-equivalent age and metallicity may be affected 
by the existence of the young component if the frosting 
model applies to elliptical galaxies. 

Finally, we make our sample publicly available to the 
community 
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TABLE 1 
Sample Definition 



Sample Size 



Description 



Lowz 
Environ 
PhotoRS 
SpecRS 
Elliptical 
Bonus 



87, 623 
57,885 
37, 026 
32, 726 
1923 
430 



All galaxies with SDSS imaging below z < 0.05 
Galaxies with Mj. < —19.0 in Lowz 
Galaxies that pass red sequence cuts (§2.1) 
Galaxies with SDSS spectroscopy in PhotoRS-Sam. 22, 621 with cr > 70 km 
The elliptical galaxy sample 
The bonus elliptical galaxy sample without SDSS spectroscopy 
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TABLE 2 

lNDEX-(T RELATION 



Index 


CoeflFecient 


All 


Rich group 


Poor group 


Field 


H/3 (A) 
H/3 (A) 


Cl 
C2 


2.02 ± 0.01 
-0.71 ± 0.04 


1.97 ±0.01 
-0.61 ±0.06 


2.02 ± 0.01 
-0.61 ±0.07 


2.08 ±0.01 
-0.81 ±0.10 


H5a (A) 
H<5a (A) 


Cl 
C2 


-1.18 ±0.01 
-4.48 ± 0.06 


-1.31 ±0.02 
-3.83 ±0.10 


-1.17 ±0.01 
-4.70 ± 0.10 


-1.03 ±0.02 
-4.79 ±0.15 


Hd> (A) 
H5f (A) 


Cl 
C2 


0.60 ±0.01 
-1.78 ± 0.04 


0.53 ±0.01 
-1.46 ±0.07 


0.61 ±0.01 
-1.88 ±0.07 


0.67 ±0.01 
-1.87 ±0.10 


H7A (A) 
H7A (A) 


Cl 

C2 


-5.08 ±0.01 
-4.67 ±0.06 


-5.20 ±0.02 
-3.99 ±0.10 


-5.09 ±0.01 
-4.80 ±0.10 


-4.99 ± 0.02 
-5.15 ±0.15 


H7F (A) 
Htf (A) 


Cl 
C2 


-1.11 ±0.01 
-2.62 ±0.04 


-1.18 ±0.01 
-2.32 ± 0.06 


-1.11 ±0.01 
-2.67 ±0.06 


-1.05 ±0.01 
-2.69 ± 0.09 


(Fc) (A) 
(Fc) (A) 


Cl 
C2 


2.59 ±0.01 
0.93 ±0.08 


2.60 ±0.01 
0.84 ±0.07 


2.58 ±0.01 
0.97 ±0.08 


2.59 ± 0.01 
1.01 ± 0.11 


Mgb (A) 


Cl 
C2 


rt.L/W _1_ L/.W-L 

3.44 ± 0.07 


3.31 ±0.07 


Q qq + n ni 
3.55 ±0.08 


q qi + n m 
3.31 ±0.11 


[MgFc]' (A) 
[MgFc]' (A) 


Cl 
C2 


3.25 ±0.01 
1.96 ±0.05 


3.28 ±0.01 
1.85 ±0.05 


3.25 ±0.01 
2.03 ±0.06 


3.21 ± 0.01 
1.97 ±0.08 


C24668 (A) 
C24668 (A) 


Cl 
C2 


6.28 ±0.02 
5.78 ±0.15 


6.33 ± 0.03 
5.51 ±0.15 


6.21 ±0.02 
6.24 ±0.16 


6.23 ± 0.03 
5.72 ± 0.23 


Ca4227 (A) 
Ca4227 (A) 


Cl 
C2 


0.99 ±0.01 
0.58 ±0.06 


1.00 ±0.01 
0.52 ± 0.05 


0.98 ±0.01 
0.63 ±0.06 


0.99 ± 0.01 
0.46 ± 0.08 



Note. — Index = ci + C2 (log^Q a — 2.2), where cr is in km s 
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TABLE 3 

SSP PARAMETBR-cr RELATION 



SSP parameters 


Coeffecient 


All 


Rich group 


Poor group 


Field 


logio Age (Gyr) 
logio Age (Gyr) 


Sl 
S2 


0.78 ±0.01 
0.50 ±0.04 


0.81 ±0.01 
0.45 ±0.05 


0.77 ± 0.01 
0.50 ±0.06 


0.74 ±0.01 
0.45 ±0.08 


[Fe/H] 
[Fe/H] 


Sl 
S2 


-0.06 ± 0.01 
0.37 ± 0.03 


-0.06 ± 0.01 
0.36 ± 0.04 


-0.06 ± 0.01 
0.42 ± 0.05 


-0.05 ± 0.01 
0.36 ± 0.06 


[Mg/Fe] 
[Mg/Fe] 


Sl 
S2 


0.25 ±0.01 
0.32 ±0.02 


0.25 ±0.01 
0.29 ± 0.03 


0.24 ±0.01 
0.35 ±0.04 


0.24 ±0.01 
0.27 ±0.05 


Note. — SSP parameters = si 


+ S2 (logio <^ - 


2.2), where a is 


in Ism s~^. 
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Fig. 1. — Randomly selected examples in the final Elliptical sample. We show the images, combined from images in gri bands, of 25 
galaxies in ascending order of velocity dispersion (a) from 70 km to 325 km s~^, from top to bottom and from left to right. 
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Fig. 2. — Examples of bulge-dominated galaxies that pass parameter cuts in preselection but fail visual inspection. Prom left to right, 
we show an SBO galaxy with a broad bar in the left panels, an SO galaxy with a faint dust lane in the middle panels and an SO galaxy 
with faded spiral arms in the right panels. Top panels: images, comlsined from images in gri bands. Middle panels: deblended images in 
r band. The yellow contours are the isophotes, and the red lines are the ellipses determined by the Ellipse algorithm. Bottom panels: 
surface brightness profile in r band, decomposed to a Sersic component (the bulge, B) and an exponential component (the disk, D). a 
stands for the major axis, e represents the ellipticity, and fir indicates the surface brightness in r band. 
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Fig. 3. — Broad-band properties of elliptical galaxies, compared with parent samples. Red points indicate the galaxies in the final 
Elliptical sample. The gray scales show the density distribution of the parent samples (LowZ in the top left panel, PhotoRS in the bottom 
left panel and SpecRS in the right panels). Top left panel: Color- magnitude diagram oi g — i vs. M^. The solid line shows the red sequence 
cut. Top right panel: Bulge-to-total (B/{B + D)) ratio vs. velocity dispersion (cr). Bottom left panel: Color gradient vs. Mr, where 
the color gradient is defined as the difference between color within 15% light radius (ris) and color between the 15% and 90% light radii. 
Bottom right panel: Color gradient vs. velocity dispersion (cr). Elliptical galaxies dominate the bright end and form a different sequence 
from the disk-dominated galaxies in the color gradient- magnitude/cr space. 
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Fig. 4. — Distribution of velocity dispersion (cr) and cllipticity (e) of red-sequence galaxies and elliptical galaxies. Top panels: Normalized 
histogram of a and e. Black solid lines indicate the distribution of the parent SpecRS sample, and red dashed lines represent that of the 
Elliptical sample. Bottom panels: Fraction of elliptical galaxies in the parent SpecRS sample, as a function of cr and e. The fraction is a 
strong function of a and e. Elliptical galaxies dominate at the massive end and at the low e end. 
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Fig. 5. — Angular distribution of the parent SpecRS sample and part of the Elliptical sample in equatorial coordinates. In the upper 
panel, we show the galaxies in SpecRS as black points. We show the elliptical galaxies in the bottom panel. Blue open triangles represent 
field galaxies, magenta filled squares represent rich-group galaxies. We also show the position of the Coma cluster in the bottom panel. 
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Fig. 6. — Distribution of velocity dispersion (tr) of elliptical galaxies as a function of environment. We show the distribution of number 
of galaxies per 20 km s~^. For clarity, we have multiplied the number of the whole sample by 50%. The field galaxies tend to have lower 
a. The median a for field, poor-group, and rich-group galaxies are 170, 186, and 203 km s~^, respectively. 
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Fig. 7. — Average spectra (unsmoothed version) of elliptical galaxies as a function of velocity dispersion (a). Before staclsing each 
spectrum, we normalize it to the mean flux between 5200 A and 5800 A where the spectrum is relatively flat. When staclting the spectra, 
we mulptiply each spectrum in each subsample by a weight factor that is the ratio of the number of the galaxies in the whole sample to 
that in the subsample at the same velocity dispersion (cr. Figure IBJ. This ensures that we are comparing spectra in all the subsamples, 
e.g., in different environments, with the same effective cr distribution. We also calculate the jackknife errors by dividing each subsample 
into 10 sub-subsamples with equal number of galaxies. The typical errors are smaller than 1% and we do not show them in this plot. We 
define the a bins to be of roughly equal size in log space. The two dotted lines are the same in each panel (0.35, 0.95) to guide the eye. 
The average spectra at lower a are apparently bluer. 
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Fig. 8. — Average spectra (uniform version) as a function of velocity dispersion (u). We show the ratio of the average spectra of all 
elliptical galaxies in each a (in km s~^) bin to that of all elliptical galaxies. 



Stellar Populations ol Ellipticals 



31 




Fig. 9. — Average spectra (unsmoothed version) of elliptical galaxies as a fu nctio n of environment. They look strikingly similar to each 
other. However, small differences do exist when we take a closer look in Figure [TOl Note it is important to control for a when comparing 
samples in different environments (See Figure [7] and |8j . 
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Fig. 10. Comparison of average spectra (smoothed version) of elliptical galaxies in different environments. We show the ratio of the 
average sjjcctra of the subsaniplcs to that of the whole sample. We calculate the jackknife errors by dividing each subsample into 10 sub- 
subsamplcs with equal number of elliptical galaxies and show them separately with gray scales. The field galaxies have a bluer continuum 
and stronger Balmcr and nebular emission lines, but by only ~ 1 percent compared to the whole sample. The rich-group galaxies, on the 
contrary, have a redder continuum and weaker Balmer and nebular emission lines, but also by only ~ 1 percent. 
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Fig. 11. — Comparison of average spectra (uniform version) of field cUiptical galaxies witii tliat of rich-group elliptical galaxies as a 
function of velocity dispersion (tr). N(ficld) indicates the number of field elliptical galaxies, and N(rich group) indicates that of rich-group 
elliptical galaxies in each a bin. The gray scales shown are jackknife errors. Compared to rich-group galaxies, field galaxies have a stronger 
continuum and stronger Balmer and nebular emission lines. The environmental dependence appears to be strong at low cr and vanishes in 
the highest a bin. 
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Fig. 12. — Fraction of galaxies with the mass fraction of the young component /y > 2.0% as a function of environment and velocity 
dispersion (a). We fit each spectrum with a two -component mode l (young+old, §4.1.4). /y is the mass fraction of the young component. 
jy = 2.0% gives NUV— r ~ 5.4, which is the cut ISchawinski et al.l ((2007) adopted to indicate recent star formation. This is based on the 
assumptions that age is the dominant factor that shapes SEDs, and elliptical galaxies consist of an old base stellar population and a small 
fraction of young stars (the frosting model). The fraction of galaxies with significant young stellar populations is a strong function of cr 
and environment, consistent with the average spectra. The fraction is higher at low a and in the field. The environmental dependence is 
strong at low cr and vanishes at high a. 
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Fig. 13. — Balmer absorption indices in the Lick indices as a function of velocity dispersion (ct) and environment. We show the 
measurements for field galaxies with blue open triangles in the left panels, those for rich-group galaxies with red upside-down triangles in 
the middle panels and those for average spectra in the right panels. We also single out the Coma galaxies as purple solid circles in the 
middle panels. The blue dashed lines represent the linear fits to the measurements for field glaxies, the red solid lines are the linear fits to 
those for rich-group galaxies, and the magenta dot-dashed lines are for poor-group galaxies (not shown here). We show the typical errors of 
measurements in the upper left corner in each panel. We correct the emission line infill by assuming EW(H/3)=0.6 EW([0 III] A5007) and 
higher-order Balmer decrement H7/H/3 = 0.46 and H(5/H/3 = 0.26. The Balmer indices are all weaker at higher a and stronger in the field. 
The emission line infill corrections of H/3 measurements of the average spectra in the lowest a bin appear to be underestimated (Figure [TSt . 
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Fig. 14. — The same as Figure [131 but for metallicity indicators in the Lick indices. The metal absorptions look very similar. Only the 
a-element indices seem to be slightly weaker in the field, but at a barely significant level. 
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Fig. 15. — Average spectra within a narrow wavelength range containing H/3, [O III], and Ha. Red solid lines represent the average 
spectra of rich-group galaxies and blue dotted lines indicate that of field galaxies. This shows the difficulty in emission line infill correction. 
When calculating the Balmer Lick indices, correcting Balmer emission lines using direct Ha or H/3 measurements merely recovers the 
template used to fit the continuum and it is in fact almost impossible to measure reliably when the emission is weak and entangled with 
absorption. Meanwhile, [O III] A5007 can be measured more reliably, but the H/3/[0 III] A5007 ratio suffers from a large scatter. 
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Fig. 16. — H/3 vs. (Fe) diagram. We plot H/3 and (Fe) on top of the S07 models with [Mg/Fe]=0.30. We here only show the galaxies 
within the velocity dispersion range 2.10 < log^o (7 < 2.30 (125 km s'^ < cr < 200 km s"!). The filled stars are the median values of field 
galaxies, poor-group galaxies and rich-group galaxies in each panel. The error bars show the Icr scatters in both axes. In the lower right 
panel, we also show the elliptical galaxies in the T05 sample, corrected to fiux-calibrated Lick indices with empirical corrections given by 
S07. Compared to rich-group galaxies, field galaxies are more spread into younger grids. The Coma galaxies in our sample appear to be 
in general older than average. Our H/3 measurements are systematically larger than those in T05 sample, most likely caused by emission 
correction. 
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Fig. 17. — H/3 vs. (Fe) diagram, but with TMB models with [Q:/Fe]=0.30. We correct all the measurements to the standard Lick indices 
with empirical corrections given by S07 to match the TMB models. In the lower right panel, we also show the elliptical galaxies in the T05 
sample. 
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Fig. 18. — H/3 vs. [MgFe]' diagram. We plot the measured H/3 and [MgFe]' on top of the S07 models with [Mg/Fe]=0.30. We only show 
the galaxies within the velocity dispersion range 2.10 < logj^Q a < 2.30 (125 km s~^ < cr < 200 km s"'^). The filled stars are the median 
values of field galaxies, poor-group galaxies and rich-group galaxies in each panel. The error bars show the la scatters in both axes. In the 
lower right panel, we also show the elliptical galaxies in the T05 sample, corrected to fiux-calibrated Lick indices with empirical corrections 
given by S07. Because the T05 catalog only includes the average index (Fe), we instead show the index [MgFe] = -^Mg b ■ (Fe). 
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Fig. 19. — H/3 vs. [MgFe]' diagram, but with TMB models with [a/Fe]=0.30. We correct all the measurements to the standard Lick 
indices with empirical corrections given by S07 to match the TMB models. In the lower right panel, we also show the elliptical galaxies in 
the T05 sample. 
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Fig. 20. — (Fe) vs. Mgfe diagram. We plot the measured (Fe) and Mgfe on top of the S07 models with age 7 Gyr. We only show the 
galaxies within the velocity dispersion range 2.10 < log^Q o- < 2.30 (125 km s'^ < a < 200 km s"^). The filled stars are the median values 
of field galaxies, poor-group galaxies and rich-group galaxies in each panel. The error bars show the l-c scatter in both axes. In the lower 
right panel, we also show the elliptical galaxies in the T05 sample, corrected to flux-calibrated Lick indices with empirical corrections given 
by S07. Rich-group galaxies are slightly more spread into higher [Mg/Fe] grids. 
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Fig. 21. — (Fe) vs. Mgfe diagram, but with TMB models with age 7 Gyr. We correct all the measurements to the standard Lick indices 
with empirical corrections given by S07 to match the TMB models. In the lower right panel, we also show the elliptical galaxies in the T05 
sample. 
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Fig. 22. — SSP-equivalent parameters vs. velocity dispersion (cr). Using H/3, (Fe) and Mg6, we derive the SSP-equivalent parameters 
from the S07 models with age, [Fe/H], [Mg/Fe] and [Z/H] = [Fe/H]+0.94 [a/Fe], The red solid lines are the linear fits to the derived 
parameters for rich-group galaxies, the magenta dot-dashed lines represent those for poor-group galaxies and the blue dashed lines are for 
field galaxies. The large open symbols show the derived parameters for the average spectra, whose derived age and [Fe/H] are affected by 
the uncertainties in emission line infill correction of H/3. The large filled symbols represent the derived parameters at the same velocity 
dispersion of the average spectra, but assuming the index-cr scaling relations (Eq. 1) hold exactly. The error bars in the upper left corner 
in the left panels are the median Monte Carlo errors. All parameters strongly correlate with cr, but all with a large scatter. More massive 
galaxies are older, more metal-rich and more strongly o— enhanced. Rich-group galaxies are systematicall older than field galaxies, by ~ 1 
Gyr, although this effect is most pronounced at low a. Rich-group galaxies also appear to be slightly more iron-poor (in terms of [Fe/H]) 
and slightly more strongly a— enhanced, but only at a barely detectable level. There is no noticeable difference of total metallicity [Z/li] 
in different environments. 



Stellar Populations of Ellipticals 



45 




Fig. 23. — C omparison of the best-fit scaUng relations between SSP-equivalent parameters and velocity dispersion (c) with those derived 
in Ncl an et al.l | |200 5) and T05. The error bars are the Icr scatter for each parameter in the whole sample, not the errors of the fits in Table 
[3] We have assumed [ci/Fe] = [Mg/Fe] and [Z/H] = [Fe/H]+0.94 [o/Fe] (TMB) in the comparison. The relations are in good agreement with 
each other. The main difference is in the age. Our derived age is ~ 3 Gyr younger than theirs, mainly because our H/3 measurements are 
larger than theirs by ~ 0.3 A (see, e.g., Figure [T6t . The difference is mostly likely caused by emission line infill correction (see text). The 
median emission correction for H/3 in our sample is 0.37 A, which translates to ~ 4 Gyr in the derived age. 



